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Abstract 

Recently, a recursion relation has been developed, generating the four-dimensional integrand of 
the amplitudes of N = 4 supersymmetric Yang-Mills theory for any number of loops and legs. In 
this paper, I provide a comparison of the prediction for the two-loop six-point maximally helicity- 
violating (MHV) integrand against the result obtained by use of the leading singularity method. 
The comparison is performed numerically for a large number of randomly selected momenta and 
in all cases finds agreement between the two results to high numerical accuracy. 



1 Introduction 



Scattering amplitudes in gauge theories are fascinating quantities partly because they provide a 
direct link between theory and experiment and partly because their investigation continues to 
uncover rich structures in quantum field theory Our understanding of computing amplitudes has 
undergone a revolution over the past decade and a half, owing in large part to the development 
of on-shell recursion relations for tree-level amplitudes [1„ 2] and a purely on-shell formalism 
for loop-level amplitudes, the modern unitarity method [3H34], These powerful modern methods 
have to a large extent made the more traditional approach of Feynman diagrams obsolete for 
tree-level and one-loop amplitudes. Thus, the current frontier is the development of systematic 
approaches for computing two-loop amplitudes. 

The maximally supersymmetric gauge theory in four dimensions, the N = 4 super Yang- 
Mills theory, has attracted a great deal of effort over the past several years, serving as a lab- 
oratory for testing new ideas. This has lead to the discovery of a surprising new symmetry in 
the planar limit, the so-called dual conformal symmetry, which, although not inherited from 
the theory's Lagrangian in any obvious way, is nonetheless a property of its amplitudes. This 
symmetry combines with the standard superconformal symmetry, dictated by the Lagrangian, 
to form an infinite-dimensional symmetry algebra, the Yangian of the superconformal group. 
This integrable structure has been intensively studied |35H50| and has recently been exploited 
successfully to determine the theory's scattering amplitudes recursively in the number of loops 
|51j . Another spectacular recent advance, following earlier insights based on a reformulation of 
the S matrix of M = 4 SYM as a contour integral on a Grassmannian manifold [52-59] — but 
this time extendable to any planar supersymmetric gauge theory — has been the development 
of recursion relations a la BCFW for the (strictly four-dimensional) loop integrand |60H63j . 

Ensuring the soundness of new ideas requires performing careful tests of the results they pro- 
duce. In this spirit, the question we wish to address in this paper is whether the result in refs. 
|6CH [63] for the two-loop six-point MHV integrand of N = 4 SYM theory can be reproduced by 
more traditional methods. In the generalized-unitarity based approach, an L-loop amplitude is 
expressed as a linear combination of known basis integrals, plus terms that are rational functions 
of the external momenta (and which will not be discussed in this paper), 

Amplitude = coefficient j Integral^ + Rational . (1-1) 

j'SBasis 

The coefficients of the integrals (evaluated in dimensional regularization) are obtained by chang- 
ing the integration range from (M 4_2e )® L into specific contours a of real dimension 4L, embedded 
in C 4i and encircling the points where the maximum number of denominators of the integrand 
become zero. Unlike the path followed in the leading singularity method |70H74] in which one 
allows any choice of contour a (by virtue of first having expanded the amplitude in an artfully 
chosen, typically overcomplete, basis), in generalized unitarity the contours a are subject to the 
constraint that any function which integrates to zero on (M 4_2e )® i must also integrate to zero 
on <7 [34] . As argued in this paper, multidimensional contours a satisfying this consistency condi- 
tion are guaranteed to produce correct results for scattering amplitudes in any gauge theory, not 
only M = 4 SYM theory. The change of integration contour has the effect of transforming the 



integrals in eq. (1.1) into contour integrals (which are easily evaluated by taking residues) . By 
making the various allowed choices of contours a one produces a set of linear equations satisfied 
by the integral coefficients which can then be solved to determine them uniquely. 
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The set of linear relations between two-loop integrals includes the set of all integration-by- 
parts (IBP) identities |75fl93j between the various tensor integrals arising from the Feynman rules 
of gauge theory; however, at present a complete knowledge of such relations involving six-point 
two-loop integrals is not available [94]. For this reason, in this paper we will carry out the analysis 
following the leading singularity method, allowing any contour a. This approach was shown to 
reproduce the result of a unitarity-based calculation |95] for the parity-even part of the two-loop 
six-point MHV amplitude in ref. [73]. By expressing the full two-loop amplitude (i.e., including 
both parity-even and odd parts) in terms of the same basis as used in ref. [73], we therefore 
expect the leading singularity method to also produce the correct result for the parity-odd part 
of this amplitude. 

The results of refs. |60[ 163] are expressed in terms of a different basis from that of ref. |73j . 
and thus it is not meaningful to check agreement between individual integral coefficients in 
either representation. A quantity that can be meaningfully compared is the two-loop integrand: 
in general, in the planar limit of any field theory, the loop integrand is a well-defined rational 
function of the external momenta (which for example can be thought of as being produced by 
the Feynman rules). We have evaluated the integrand for a large number of randomly selected 
momenta and in all cases find agreement with the recent literature to high numerical accuracy 

An interesting spinoff of the calculation in this paper are the potential applications of the 
intermediate results (in particular, the enumeration of the global poles of the integrand and the 
expressions for the cut integrals). Once all the necessary IBP relations do become available, we 
expect these results to greatly facilitate the task of determining the maximal-cut contours that 
allow the extraction of integral coefficients in any two-loop six-point gauge theory amplitude. 

A complementary approach is that of refs. [64H68] . in which the heptacut integrand is recon- 
structed by polynomial matching in similarity with the OPP approach |21j . A recent paper by 
Zhang |69] adds tools required in such an approach for reducing integrands to a basis of mono- 
mials. 

This paper is organized as follows. In Section [2] we explain conventions and introduce nota- 
tion used throughout the paper. In Section [3] we compute the heptacuts of the general double 
box integral and of the factorized double box with six massless legs. In Section [4] we apply these 
heptacuts to the ansatz for the two-loop six-point MHV amplitude of N = 4 SYM and employ 
the leading singularity method to obtain linear equations satisfied by the integral coefficients. We 
discuss how these equations can be used to directly obtain the parity-even part of the integrand. 
We then explain how to evaluate numerically the full N = 4 SYM integrand (i.e., including both 
parity-even and odd parts) as obtained by the leading singularity method and report agreement 
with the result in refs. |60[ 163] , Finally, in Section [5] we provide conclusions and suggest direc- 
tions for future investigation. In Appendix [A] we provide full details of all the heptacuts of the 
two-loop six-point MHV amplitude. 

2 Notation and conventions 

In this section we explain conventions and introduce notation used throughout the paper. 

All external momenta in an amplitude are outgoing and will be denoted by ki. We will 
make use of the spinor helicity formalism |97H104| in which a given massless four-dimensional 
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momentum is written as a tensor product of two massless Weyl spinors, 
We define the spinors 

Xi = u + (ki), Xi = u-(ki) 
and the Lorentz invariant inner products formed out of the spinors, 



(ij) = (i \j > = u-(ki)u + {kj) 



[ij] = (i \j > = u + (ki)u-(kj) 



which satisfy 



Spinor strings are defined as follows 

(Kr\f + Q\Kj 



(K~ 



(Kr\f\KJ) + (KT\Q\KJ) 
(KiP^PKj] if P 2 = 0. 



(2.1) 
(2.2) 

(2.3) 
(2.4) 

(2.5) 
(2.6) 



Si 



We will use the following notation for sums and invariant masses of external momenta, 

= fcjj + • • • + ki n (2-7) 

(kn + • ■ ■ + h n f (2.8) 

Kf . (2.9) 

Throughout we will make use of the "flattened" momenta introduced in ref . |21|. [23] : for a pair 
of momenta K\, K 2 , define the quantity 

11 ± = {K 1 -K 2 )±^, Ai = (Kx ■ K 2 f - K\K\ (2.10) 

which can take two different values if both momenta are massive (i.e., if S\S 2 7^ 0). For a given 
value of 71 one defines a pair of massless "flattened" momenta as follows 



K 



Ki - (Si/7i±)#2 



i± 



1 



SiS 2 h 2 



K 



2± 



1± 



K 2 - {S 2 / 11± )K 1 
1 - 5i5 2 / 7l 2 ± 



(2.11) 



If one of the momenta K\ or K 2 is massless, 7i± can only take one value, and we will use the 
following abbreviated notation: 



S\S 2 
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Similarly, we will use the notation 

72 ± = (K 4 ■ K 5 ) ± 



A 4± 



K 4 - (S 4 / 72 ±)K 5 
1-S 4 S 5 /jI ± 



A, 



K\ 



72 



S4S5 



K 
K 



2K X ■ K 2 

K t - {Sxhx)K 2 

K 2 - (S 2 h l )K l . 

= (K 4 ■ K 5 ) 2 - K\Kl 
K 5 - [S b h 2 ±)K A 
5± ~ l-S A S,h 2 ± 

2K 4 ■ K 5 

K 4 - (£4/72)^5 
K 5 - (S 5 / 12 )K 4 . 



(2.12) 

(2.13) 
(2.14) 

(2.15) 



Finally, we will denote the elements of the dihedral group Dq as follows 

(H = (1,2,3,4,5,6) a 2 = (2,3,4,5,6,1) * 3 = (3,4, 5,6, 1,2) 

a 4 = (4,5,6,1,2,3) a 5 = (5,6,1,2,3,4) ^ = (6,1,2,3,4,5) 

a 7 = (6,5,4,3,2,1) a 8 = (5, 4, 3, 2, 1, 6) ^ = (4,3,2,1,6,5) 

(710 = (3,2,1,6,5,4) a n = (2,1,6,5,4,3) a 12 = (1, 6, 5, 4, 3, 2) 



(2.16) 



3 



3 Heptacut two-loop integrals 



As preparation for computing the generalized-unitarity cuts of the two-loop six-point amplitude, 
in this section we compute the heptacuts of the general double box integral and the factorized 
double box integral, illustrated below in Figs. [T] and [5] These integrals are part of the linear bases 
in which the two-loop six-point N = 4 SYM helicity amplitudes were expanded in ref. |73| 195) 
1105] . In this paper we shall use the (over complete) basis in ref. [73j (illustrated for convenience in 
Fig. [7]in Section]!]). As explained further in Section 4.1 the heptacuts of the remaining integrals 
in this basis are easily obtained from the heptacuts of the double box integrals by multiplying 
additional factors (involving propagators and possibly numerator insertions) onto the integrand 
of the latter. From the knowledge of the generalized cuts of all basis integrals, the cuts of the 
amplitude are easily obtained, as will be explained further in Section [4j 

We emphasize that the generalized cuts considered throughout this paper are strictly four- 
dimensional (as opposed to (4 — 2e)-dimensional). Moreover, we will suppress the Feynman ie- 
prescription when writing propagators. Finally, the internal lines in all diagrams are taken to be 
massless. 



3.1 The maximal cut of the general double box 

In this section we will compute the maximal cut of the double box integral 




Figure 1: The general double-box integral. 



where the vertex momenta K\ shown in Fig. [T] are typically sums of several lightlike momenta, 
Ki = ki x + ■ ■ ■ + ki n with /c?. = 0. Setting D = 4 — 2e, this integral is defined by 

d D h d D £ 2 f 1 1 1 1 1 1 1 

(2vrp (2xp K x f (h -K x - K 2 f {l x + £ 2 + Ktf Q (£ 2 - K h f (£ 2 -K 4 - Ktf 

(3.1) 

where the integration is over real Minkowski space for both loop momenta. By considering the 
heptacut of the general double box integral we can easily obtain the heptacuts of all the double 
boxes appearing in the basis in Fig. [7] (in Section [1]) by taking appropriate vertex momenta to 



be massless. As will be explained in Section 3.1.5, in the case of six massless external momenta, 



there are three qualitatively distinct ways of distributing the momenta at the vertices of the 
double box. In order to streamline the presentation we will first compute the heptacut of the 
completely general double box integral, making no assumptions about masslessness of any vertex 



4 



momenta. Subsequently, we will describe each of the three cases in turn. 



Formally speaking, the four-dimensional heptacut of (3.1) is obtained by replacing each of the 
seven propagators by a 5-function whose argument is the denominator of the propagator in ques- 
tion, and replacing the D-dimensional integration measure by the corresponding four-dimensional 
measure (up to factors of 2ir, depending on conventions). Thus, the heptacut replaces the double 



box integral in (3.1 ) by the integral 

Jformai = J dH x d% 5{£ 2 ) &((l x - K^ 2 ) 8[{h - K x - K 2 ) 2 ) 5((£ l + £ 2 + K 6 ) 2 ) 

x 8{i 2 2 ) 5{(£ 2 - K 5 ) 2 ) 5((£ 2 -K A - K 5 ) 2 ) . (3.2) 

However, this integral only receives contributions from regions of integration space where the 
loop momenta solve the joint on-shell constraints 





= 


(3.3) 


(h - K t ) 2 


= 


(3.4) 


(h -Kx- K 2 f 


= 


(3.5) 


a 


= 


(3.6) 


{h - K 5 ) 2 


= 


(3.7) 


{h ~K A - K 5 ) 2 


= 


(3.8) 


{h + £ 2 + K e ) 2 


= o, 


(3.9) 



which in general only have solutions for complex loop momenta (£ x ,£ 2 ) 6 C 4 x C 4 

The natural definition of ^-functions with complex arguments involves contour integrals - 
integrating out a variable q in an integrand involving <5-functions will fix q to some value qo; i n 
the language of contour integrals, this corresponds to integrating in the complex g-plane along 
a small circle centered at q^. Indeed, as observed in ref. [106\ 1107] . Cauchy's residue theorem 
implies that the localization property 

dqS(q-q )f(q) = f(q ) (3.10) 

remains to hold if we define S(q — qo) = ^^z^ and take the integral to be a contour integral 
along a small circle in the complex g-plane centered at go- 

By analogy, taking the four-dimensional heptacut of the double box integral ( |3.1[ ) should really 
be understood as a change of integration range from IR^ x M. D to a surface (of real dimension 8) 



embedded in C 4 x C 4 while leaving the integrand in eq. (3.1 ) unchanged. The maximal-cut integral 
is thus a multidimensional contour integral whose contour is in general a linear combination of tori 
encircling the so-called global poles of the integrand. These are points (£i,£ 2 ) 6 C 4 x C 4 where all 



seven propagators in (3.1) become on-shell Q The change of contour away from real Minkowski 



space renders the double box integral IR and UV finite, and one can therefore disregard the 
dimensional regulator part of the measure d~ 2e £id~ 2e £ 2 and the (— 2e)-dimensional components 
of the loop momenta. 

In the following we will continue to write multidimensional contours symbolically in terms of 



5-functions, as in eq. (3.2), as we find the latter notation more suggestive. As it turns out, in all 



1 To be more accurate, we will use the terminology "global poles" to refer to points where, in addition to the 



seven on-shell constraints ( 3.3 H 3.9 1 being solved, there is an additional singularity, for example coming from the 



Jacobians arising from linearizing these cut constraints. A more precise definition will be given in Section 4.1.3 



5 



cases considered in this paper, the only respect in which the multidimensional contour integrals do 
not behave like integrals of 5-functions is the transformation formula for changing variables: Given 
a holomorphic function / = (/l, • • • , /n) : C n — > C n with an isolated zercj^jat a £ C n , the residue 
at a is computed by the integral over the contour T e (a) = {z £ C n : = Ci, i = 1, • • • , Ji}. 

This contour integral satisfies the transformation formula 

1 f h(z)dz 1 A • ■ • A dz n _ h(a) 
(2*i)» ./r.(.) h{z)---f n {z) - detail 1 ' ' 

which, crucially, does not involve taking the absolute value of the inverse Jacobian. This ensures 
that this factor is analytic in any variables on which it depends, so that further contour integra- 
tions can be carried out. 



In order to visualize the multidimensional tori in question, it turns out to be convenient to 
use the following parametrization of the loop momenta: 



ocxKY + a 2 K^ + a 3 (K\-\^\K b 2 
+ {3 3 (Kt\^\K b 5 - 



) + a 4 (K b 2 -\^\K[ 



l^K 



(3.12) 
(3.13) 



The virtue of this parametrization is that it linearizes as many of the cut constraints as possible, 
and in turn it becomes easy to locate the positions of the global poles of the integrand in 
the coordinates ai,..., a 4 , Pi, ... , j3&. The multidimensional tori discussed above are then easily 
obtained as products of small circles each encircling one of the entries of a given global pole. 

After changing variables from the components of the loop momenta £ x and l\ to the param- 
eters o-i and f3j, the heptacut of the double box integral becomes 
i 



J 



1 



7i7| 



=i 

x 5 

x 5 
x5 



|U,i OOLi 



, d£ 2 

v, 3 d(3j 



5(aiCt2 



((A 



71 

12 



«2 



ill 

71 



40304) dy(ai — 1) ( 02 
4a 3 a 4 ) SiPifc 



Si 
71 



40304 



4/3 3 /3 4 ) 



(02 



4/3 3 /?4 )8 [{/3 



12 



1 



Ih 



S4 

72 



4/3 3 A 



+ £2 + K 6 f 



Iparam . 



(3.14) 



where the subscript "param" on the argument of the 5-function in the last line indicates that it 
is to be evaluated in the parametrization ( |3.12 )-(3.13) and where the Jacobians associated with 
the change of variables are, respective!}^] 



det * 

AM OOLi 



ill 



det id 

u,j dfij 



ill 



(3.15) 



In the parametrization (3.12)-(3.13), the six on-shell constraints (3.3)-(3.8) which only involve a 
single loop momentum are solved by setting 



TiQSj + 7i) 

ll-S!S 2 ' 

545*5(55 + 72) 
72 (5 4 5 5 -7I) 



OL2 



ft 



5i5 2 (5i + 71) 

7i(Si5 2 -7i 2 ) 
_ 72(5 4 + 72) 



ll 



OLZOL^ 



/3 3 /3 4 



5i5 2 (5i + 7i)(5 2 + 71) 

4( 7l 2 - 5i5 2 ) 2 
5 4 5 5 (5 4 + 72) (5 5 + 72) 



4(7l 



(3.16) 
(3.17) 



5 4 5 5 ) 2 

iff by choosing a small 



2 A function / = ...,/„) : C" — > C" is said to have an isolated zero at a G C 
enough neighborhood U of a one can achieve / _1 (0) PI U = {a}. 

3 These identities can be shown by writing the determinants as square roots of Gram determinants (since both 
of these equal the volume of the spanned parallelotope, up to a complex phase) and using the special properties of 



the vectors {K^, K^ 1 , {K[~\~/' 



\K b 2 -),(K>- 



\r 



\Kt)} and {K b f,K b ^,(K 



^\Kf),(KZ-\r\KT)}. This will 



fix the determinants up to an overall factor of i which can be found numerically. 
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We observe that the integrations of the 5- functions in (3.14) unambiguously fix the values of 
a%, «2, Pi, p2- After imposing the last on-shell constraint (3.9), the four variables as, a.4, /?3, /?4 
are subject to three relations, and one is free to choose either one of them to be an unconstrained 
free complex parameter z. Each of these four choices give rise to different contributions to the 
heptacut of the double box. As will be explained further in Section 3.1.5, these contributions 
correspond to the various existing classes of solutions to the on-shell constraints (3.3)-(3.9), and 
the heptacut (3.14) is an appropriately weighted sum of these contributions. 

In Section |3.1.1 below we compute in detail the contribution to the heptacut obtained by 
letting z = «3 be the unconstrained parameter. The results for the remaining contributions 
are quoted in Sections 3.1.2f3~T~4 and are obtained in an entirely analogous way. Finally, in 
Section 3.1.5 we explain how to assemble the heptacut of the double box integral from its various 
contributions. 



3.1.1 Leaving z = 03 as the free parameter 



In this example, we aim to leave z = 03 as the unconstrained parameter, and we will therefore 
integrate out ol\, a-i, and P2, P3, fi/i. This will proceed in three stages: first we integrate out 
the three 5-functions involving only the a- variables; then we integrate out the three (^-functions 
involving only the /3-variables; finally, we integrate out the remaining /3-variable. 

Thus, we start by considering the following integral whose integrand consists of all the 5- 



functions in (3.14) that only involve the a- variables, 



J a = J doi\daida/± b(a.\OL2 — 40304) <5^(ai — 1) ^02 — — 40304 

x <$((ai - f - l) (a 2 - f - l) - 4a 3 a 4 ) . (3.18) 



This integral is the inverse Jacobian associated with the change of integration variables from the 
o-parameters to the arguments of the <5-functions. It is straightforwardly evaluated to yield 



J 1 a 



1 



4 1 



(3.19) 



where we recall that the (^-functions in eqs. (3.14) and (3.18) are a short-hand notation for multi- 



dimensional contour integrations. As these are subject to the transformation formula in (3.11), 



which involves the Jacobian of the transformation rather than its absolute value, there is no 



absolute value in eq. (3.19). 



Second, we consider the following integral whose integrand consists of all the 5-functions in 



(3.14) that only involve the /3-variables, 



72 



1 A 



54 
72 



4/3 3 A 



1 



(3.20) 



In addition to integrating out /3\, 02, we have here chosen to integrate out ^4. Alternatively, we 
could have imagined integrating out ^3; but this would ultimately lead to the same final result. 
Once again, the integral Jg is a Jacobian and is straightforwardly evaluated to yield 



J, 



4 ( 1 - ^# ) & 



(3.21) 
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Putting together the partial results in eqs. (3.15), (3.19), (|3.21 ) and applying the loop momentum 

I par am 



parametrization to the argument of the remaining factor 8( (li + £ 2 + Kq) 2 \ ) in eq. (3.14) 
one finds the expression 



7172 



where 
B x 

B 

B-i 



32( 7l 2 -5 1 5 2 )(7 2 2 -S 4 55)y a 3 /3 3 
(Kl-\ llx \Kl-)( ai K^ + a 2 K^ + a 3 (K\-\^\K' 2 -} + a 4 (K^\K\-) + Kg 



x (a 1 ^ + a 2 ^ + a 3 (^l7 M l^) + «4(^|7 M l^) + ^6) " \ S $ 

s 4 s 5 (s 4 + l2 ){s b + 72 )(kI-\^\kI) 

4( 72 2 - S 4 S 5 ) 2 

x (axK^ + a 2 K^ + a 3 (^~|7 M |^2~) + a 4 (K^ \^\K\-) + . 



(3.22) 
(3.23) 
(3.24) 

(3.25) 



We can integrate out /3 3 using 
In conclusion, leaving z = a 3 unconstrained and integrating out the remaining loop momentum 



6{B^ 3 + B + B^p^) 



{Bl-AB.B^Y 112 



(3.26) 



parameters from (3.14) produces the following contribution to the heptacut of the double box, 

-1/2 



J 



7172 



j> d ^[B Q {zf-AB 1 {z)B^ 1 {z) 



(3.27) 



'*=«* 32( 7l 2 -5 1 5 2 )( 72 2 -5 4 S 5 ) 
where we have relabeled a 3 = z and made the dependence of Bi,Bq, B_\ on z explicit. Although 
the factor ( • • • ) ^ 2 appears to have a branch cut, the radicand turns out to be a perfect square in 



all cases considered in this paper, and the integrand in eq. (3.27) always takes the form 



1 



ziz-P)- 



As we will discuss further in Section |4j we will allow the integration contour in eq. (3.27) to 
encircle any individual singularity of the integrand. 



Note that the form of the final result (3.27) does not depend on the order of integration: if one 



instead chooses to integrate out a\, a> 2 , 04 and fix, /3 2 , /3 3 from (3.14) and subsequently integrates 



out /?4, one finds (3.27) as may easily be checked. 



3.1.2 Leaving z = a 4 as the free parameter 



Integrating out a\,a 2 ,a 3 and f}\, /3 2 , /? 3 , f3 4 from (|3.14|) produces 

7i72 



J 



32( 7l 2 - 5!5 2 )( 72 2 - S 4 S 5 



j^B^zf-ABl^BUiz) 



-1/2 



(3.28) 



where 

Bl(z) = {Kl-\^\Kl-)( ai KY + a 2 K^ + a z {K^^ (3.29) 
B'(z) 



hK^ + fcK^ + K^ 

x (axK^ + a 2 K^ + a 3 {K\-\-f\K%-) + z{k\-\^\k[-) + K^j - (3.30) 
S 4 S 5 (S 4 + l2 )(S b + 72)^5" \l»\K\~) 

4( 72 2 - s 4 s b y 

x (axK^ + a 2 K^ + a 3 (K^\^\K^) + z(K^\^\K^) + K%) . (3.31) 
Again this result is independent of the order of the integrations. 



8 



3.1.3 Leaving z = fa as the free parameter 



Integrating out a\, a 2 , 03, 04 and j3\, (3 2 , 04 from (3.14) produces 

7i 72 



J 



32( 7l 2 -5 1 S 2 )(7 2 2 -5 4 55) 



! ' :/ ioW 2 -4iiWi-iW) V2 (3.32) 



where 



= <^-| T ^|^-) + /5 2 ^ + ^<^-| Tm |^-) + /? 4 (^-| Tai |^-> + i^ 6/ ,) (3.33) 



A (z) 



A-i(z) 



aiK^ + a 2 K^ + 

x + f3 2 K b 5fM + z(^"| 7At |^-) + fa{Kl-\ ltl \K\-) + tf 6/1 ) - §S 6 (3.34) 

S\S 2 {Si + 7i)(5 2 + 7 i)(K2~|7 M l^~) 
4( 7l 2 - 5i5 2 ) 2 



x (fti^ + 0^ + | 7M |^r) + HKt\lAKl~) +K^). (3.35) 

Again this result is independent of the order of the integrations. 



3.1.4 Leaving z = ^4 as the free parameter 



Integrating out 01,0^03,04 and /3i,/3 2 ,/?3 from (3.14) produces 

71 72 f dz ( sQ ..„,,.„ , A" 1 / 2 



J 



z =/ 3 4 



32(7 2 -5iS 2 )(7 2 -5455) 



j^AKzf-AA'^AUiz] 



(3.36) 



where 



Ai(z) = (^-|7^|^-)(^^ + /3 2 ^+ .i ! (A-;-| v |^-) + -.(/vrivlA'D + AV) (3.37) 



+ o 2 i^ + 

x (ft^ + /3 2 i^ + fo{K\-\ ltl \Ki-) + z{K\-\^\K\-) + K^) - \Ss (3.38) 
5i S 2 (S l + 71 ) (5 2 + 71 ) < ^ 2 ~ I r I ) 



4(7i 2 -SiS 2 ) 2 

x [PtK^ + /3 2 i^ + ^{K\-\^\K\-) + z(Kl~\^\Kl) + K 6fl 



(3.39) 



Again this result is independent of the order of the integrations. 



In the section below we will explain how to assemble the heptacut double box integral from 



the individual contributions in eqs. (3.27), (3.28), (3.32), (3.36). 



3.1.5 Assembling the heptacut double box from its contributions 

As alluded to in the beginning of this section, there are three qualitatively distinct ways of 
distributing six massless external momenta at the vertices of the double box. In this section 
we shall describe the classification of these cases, and how to assemble the heptacut double 



box integral from the individual contributions in eqs. (3.27), (3.28), (3.32), (3.36) in each of 



these cases. Before proceeding to describe the classification of the three cases, it is useful to 
appreciate that the individual contributions to the heptacut arise from the different existing 



classes of solutions to the joint on-shell constraints (3.3)-(3.9). Indeed, these seven constraints 
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leave one unfixed degree of freedom z in the two loop momenta £i,£ 2 , and in general there are 
several classes of solutions. Each class is parametrized by a free complex variable z and has the 
remaining seven loop momentum parameters at, f3j fixed to specific values. Two classes S and S' 
are identical if and only if there exists an invertible holomorphic function <p(z) that maps one 

class into the other; that is, ip = a i\gi anci V? i,Pi\s) = Pi\s' ^ or * = ' ' ' 

The classification into the three cases is given as stated below. 

• Case I: all three vertical propagators in Fig. [T] are part of some three-particle vertex. There 



are six classes of kinematical solutions to the on-shell constraints (3.3)-(3.9), illustrated in 
Fig. [2] The heptacut double box integral is J = X^=i Ji wnere 





Jl 


= J 


z=/3 3 


■h 


= J e 


= J 


z=a 3 




■h 


= J 




Jl 


= -h 


= J 


Z=Ct4 



(given in eq. (3.32)) 



(given in eq. (3.27)) 



(given in eq. (3.36)) 



(3.40) 



(given in eq. (3.28)) . 



This case encompasses heptacuts #1-5 and #8, discussed in detail below in Sections |4.1| 
A. 2 A. 5 and A.8| where one can also find explicit results for the on-shell values P\,Q\ 
etc. quoted in Fig. [2] To determine, for example, the function /3s(z) in solution 1S5, one 
expresses the loop momenta appearing in the on-shell constraint \[£\ + £2 + Kq) 2 = in 
their parametrized form (3.12 )-(3.13 ), sets the values of the parameters c^-, /3j equal to those 



quoted in eqs. (3.16)-(3.17) and below solution 1S5 in Fig. [2] and then solves the on-shell 
constraint for (3%. 

Case II: the left- and rightmost vertical propagators in Fig.[T]are part of some three-particle 
vertex, but the middle one is not. There are four classes of kinematical solutions to the 
on-shell constraints (3.3)-(3.9), illustrated in Fig. [3] The heptacut double box integral is 

J = YU=i J* wnere 



Jl 


= J 


z=/3 3 


a 


riven in eq. 


(3.32)) 


■h 


= J 


z=f3 4 


a 


;iven in eq. 


(3.36)) 


■h 


= J 


Z=Ct4 


a 


aven in eq. 


(3.28)) 


Jl 


= J 


z=a 3 


a 


;iven in eq. 


(3.27)) 



(3.41) 



This case encompasses heptacut #6, discussed in detail below in Section |A.6| where one 
can also find explicit results for the on-shell values Pi, Q\ etc. quoted in Fig. [3] 

Case III: the two rightmost vertical propagators in Fig. [T] are part of some three-particle 
vertex, but the leftmost is not. There are four classes of kinematical solutions to the on- 
shell constraints (3.3)-(3.9), illustrated in Fig. [1] The heptacut double box integral is J = 
Yji=i Ji where 

Jl = J 2 = J\ z=p3 (given in eq. fl3.32| )) 
J3 = J4 = J\ z=0i (given in eq. ( |3.36| >) . 
This case encompasses heptacut #7, discussed in detail below in Section A. 7| where one can 
also find explicit results for the on-shell values P^, etc. quoted in Fig. [4] In this case, 
because both K\ and K2 are massive, there are two solutions for 7^ (given in eq. (2.10)) 
and therefore two pairs of flattened momenta (K\±, K\^) (see eq. (2.11)). The on-shell 
values of the loop momenta are independent of which sign is chosen, 



(3.42) 



£i(a^ , a~2i Qjjj , cej) 
Wi + ,/3 2 + ,/3 3 + ,/? 4 + ) 



(3.43) 
(3.44) 
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Solution Si, obtained by setting 
as = P* h = z 




Solution ^3, obtained by setting 
a 3 = p 3 = 

«4 = Pi 04 = Z 

K 3 




Solution £5, obtained by setting 
as = 3 = 3 (z) 



«4 



01 








Solution 52, obtained by setting 

"3 = z 03 = <2* 

04 = 04 = 




Solution 54, obtained by setting 

a 3 = 3 = 
a 4 = z 4 = Qi 

^3 




Solution 5c obtained by setting 



a 3 = z 
04 = 



/3 3 = 
04 = 



Figure 2: The six kinematical solutions to the heptacut constraints for the double box topology 
in case I. For all solutions, the loop momentum parameters (a±, at2, 0i, 02) are set equal to the 



values given in eqs. (3.16)-(3.17). Any blob connecting more than three legs does not have a 



well-defined chirality and its sign should be ignored. For 5s and Sq, the parameters 03 and 04 
are determined by solving the on-shell constraint {l\ +£2+Kq) 2 = for the respective parameter. 
The on-shell values P\,Q\ etc. are functions of the external momenta; examples may be found 
in Sections |OJ |A~2]|A~5| [A~8 
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Solution Si, obtained by setting 

_ (fcri*6iifcrx*-Q?) a _- 

"3 - 2{k 1 k 4 )[k 5 k 2 ](z-Ql) P'i- z 
a.4 = /?4 = 




Solution 52, obtained by setting 



a 3 = 



Qi4 



2[fclfc 4 ](fc5fc2>(z-Q5 








Solution 53, obtained by setting 
„ _ n R - Q'(P2-P6)(z-Pi) 

" 3 " PZ ~ (P 2 -Pi)(z-P a ) 
«4 = Z Pi = 




Solution 54, obtained by setting 



a 3 
04 



2; 




03 

Pa 





Q2(p 2 '-p ( r)( z -p 1 *) 



Figure 3: The four kinematical solutions to the heptacut constraints for the double box topology 
in case II. For all solutions, the loop momentum parameters (ai, 0.2, Pi, P2) are set equal to the 



values given in eqs. (3.16 1-(3.17). The on-shell values P\,Q\ etc. are functions of the external 



momenta and may be found in Section A. 6 
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Solution Si, obtained by setting 



o 



± 



o 



± 



Pi 



If* 



ft = 




Solution 52, obtained by setting 



a 



Pf QJ* z-Q 



ft 
ft 








Solution 53, obtained by setting 



n 



± 



o 



± 



± 



ft 








Solution 54, obtained by setting 







a 



P-^ Q 2" 2> ■ Q 2 



ft 
ft 







Figure 4: The four kinematical solutions to the heptacut constraints for the double box topology 
in case III. In this case, because both K\ and K2 are massive, there are two solutions for 7^ 
and therefore two pairs of flattened momenta {K\ ±j K\±)- For all solutions, the loop momentum 



parameters (a\, «2, ft, ft) are set equal to the values given in eqs. (3.16 1-(3.17|. The on-shell 
values P^, Q2 etc. are functions of the external momenta and may be found in Section 
quantity ^ is defined as 



A.7 



The 



± 



5i5 2 (5 1 + 7l ± )(5 2 + 7l ± ) 



4(7? 



±2 



SiS 2 i 



(3.45) 
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3.2 Heptacut of the factorized double box 



For the double box integral considered in Section |3.1| there are various ways of distributing six 
external momenta at the vertices, and we therefore computed the heptacut with an arbitrary 
number of external legs. In contrast, the factorized double box integral 



d 



D, 



1 



f f d D £ 2 1 I 1 I 



h) 2 (h - M 2 (*2 - fc456) 2 



(3.46) 



with which we shall be concerned in this section, admits a unique way of distributing six (cycli- 
cally ordered) external momenta at its vertices, and we therefore restrict ourselves to this case. 
The factorized double box with six massless legs is illustrated in Fig. [5] below. 




fa h 

Figure 5: The factorized double box integral with six massless external momenta. 



We consider the heptacut defined by imposing the following joint on-shell constraints 





= 


(3.47) 


(h - fcl) 2 


= 


(3.48) 


(h - M 2 


= 


(3.49) 




= 


(3.50) 




= 


(3.51) 


& ~ k 56 ) 2 


= 


(3.52) 


(h - h 56 ) 2 


= . 


(3.53) 



This is not a maximal cut: we deliberately leave one degree of freedom z in the loop momenta 
unfrozen to make manifest the singularities in the Jacobians arising from changes of variables. 
This in turn makes the global poles of the integrand easy to identify. 
We will use the following parametrization of the loop momenta 

^ = aik% + a 2 kt^ + a z {k^\^\k 2 ) + a^^k^) (3.54) 
t 2 = /3i^ + /3 2 ^ + /3 3 (A: 5 -|7 M |A: 6 -)+/3 4 (A: 6 -|7lA: 5 -). (3.55) 
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The constraints ( 3.47 )-( 3.52) form a special case of eqs. (3.3)-(3.8) after appropriately relabeling 



the external momenta. From eqs. (3.16)-(3.17) we then find that the constraints are satisfied by 
setting 

a± = 1 , 02 = , 0304 = 

ft = , ft = 1 , ftft = . 



(3.56) 



The final on-shell constraint (3.53) combined with eq. (3.56) turns out to have four classes of 



solutions which can be conveniently expressed by defining the spinor ratios 



Qi 



(45) 
2(46) ' 



Qi 



[45] 
2 [46] ' 



(3.57) 



Thus, the four classes of solutions to the joint on-shell constraints (|3.47|)-(|3.53|) are 

Sr. 



5,: 



":; 


= z , 


ft 


= Qi 


Qi 


= 0, 


ft 


= 


«3 


= 0, 


ft 


= Qi 


Q'i 


= z , 


ft 


= 



5 2 : 



5 4 : 



03 


= z , 


ft 


= 


Q!4 


= 0, 


ft 


= Qi 


03 


= 0, 


ft 


= 


O4 


= z , 


ft 


= Qi 



(3.58) 



where the parameters a%, a 2 , ft, ft are put equal to the values quoted in eq. (3.56). The solutions 
are illustrated below in Fig. [6] 

Since we are not cutting all propagators, the relevant quantity to compute is not the heptacut 
of the factorized double box integral itself, but rather the Jacobian 

J = (jdH 1 5{g{)5{{i 1 -k l f)S{{l 1 -k l2 f) 

x (J d%5(£l)5((£ 2 -k 6 ) 2 )5{(h-h 6 ) 2 )5((£ 2 -k 456 ) 2 )) (3.59) 

since from this one can easily obtain the heptacut of any given object, as we will see below. The 



5-functions in eq. (3.59) may be integrated out along the lines described in Section 3.1.1 and 



one finds that for all four kinematical solutions Si the contribution to the Jacobian J is given by 

1 f dz 



J 



16S1 2 S4 5 S 56 J r . z 
Then the heptacut factorized double box integral is 

4 



for i = 1 , . . . , 4 . 



(3.60) 



— — V 

I6S12S45S56 JVi 

1=1 L 



dz 1 
V (h - k 123 ) 2 



(3.61) 



Si 



where the subscript ( • • • ) | 5 . indicates that the propagator is to be evaluated in the parametriza- 



tion (3.54)-(3.55 ) with the parameters set equal to the values in eqs. (3.56) and (3.58). Likewise, 



the heptacut (3.47)- (3.53) of the two-loop amplitude is 

1 A r dz ' 



16S12S45«56 



=1 Jt * z i=i 



(3.62) 
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Solution S±, obtained by setting 
«3 = z , h = Qi 

«4 = , /?4 = 



Solution 52, obtained by setting 



«3 = 2 
«4 = 



/3 3 = 





Solution 53, obtained by setting 
as = , pa = Q\ 



Q!4 = Z , 



/5 4 = 



Solution 54, obtained by setting 
a 3 = , ft = 



«4 = Z , 



Figure 6: The four kinematical solutions to the heptacut constraints given in eqs. ( 3.47 )-( 3.53 ) 



for the factorized double box with six massless external momenta. The on-shell values Q\ and 



Q' are defined in eq. (3.57). 
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4 Global poles of the two-loop six-point integrand 



In this section we apply the heptacuts discussed in Section [3] to compute the two-loop six-point 
MHV integrand in N = 4 SYM theory. The parity-even part of this integrand was first calculated 
in ref. [95] and was later reexamined using the leading singularity method in ref. |73j . In Section 



4.1 we give a pedagogical review of the use of the leading singularity method to determine 
the two-loop six-point integrand. In particular, we discuss how one can set up linear equations 
to determine the parity-even part of the integrand directly. The approach here is similar to 



that of ref. |73j . but differs in the use of the loop-momentum parametrization (3.12 )-(3.13) and 



(3.54 )-( 3.55 ) which has the virtue of making the multidimensional contours associated with the 
heptacuts completely explicit. It is also worth remarking on the similarity between the approach 
followed here and that of Forde in ref. [23j: the heptacut at two loops, in analogy with the triple 
cut at one loop, leaves one a priori undetermined contour integration. The associated contour may 
be chosen to encircle the various poles coming from the measure or from additional propagators, 
allowing the determination of the integral coefficients by use of the residue theorem. 

The main result of this paper is contained in Section |4.2| in which we report on a numerical 
check that the result produced by the leading singularity method for the full (i.e., parity-even and 
odd) two-loop six-point MHV integrand is in agreement with recent predictions in the literature 

[501152]. 

In general one can apply integral reductions to any two-loop six-point amplitude to express 
it as a linear combination of integrals in some sufficiently large basis. Due to the special symme- 
tries of N = 4 SYM theory amplitudes, it is natural to include integrals in the basis that reflect 
these symmetries (in particular, dual pseudoconformal integrals), and in this paper we shall use 
the basis in ref. |73j . for convenience illustrated in Fig. [7] below. Note that in order to express the 
parity-odd part of the amplitude, this basis contains additional (non-dual conformal invariant) 
integrals beyond those included in the bases of refs. |95[ 1105] and as such is over-complete. We 



shall discuss the consequences of the over-completeness in detail in Section 4.1 

Thus, we will use the following ansatz for the planar two-loop six-point MHV amplitude of 
TV = 4 SYM theory 



^6,MHV — a ^ ' (4-1) 



4 

j=l,...,24 
j=l,..,12 



in which it is expressed as a linear combination of the basis integrals illustrated in Fig. [7] whose 
symmetry factors r\ are given as follows, 



{ri,r 2 ,r 3 ,r 4 ,r 5 ,r 6 ) = (|, 1, \, \, 1, l) 

(4.2) 



(r 7 ,r 8 ,r 9 ,rio,rii,ri2) = 1,1,1, \ 



(ri3,ri4,ri5,ri6,ri7,ri 8 ) = (|, \, |, |, |, l) 

i 1 1 1 A 



(ri9,r 2 o,?"2i,r22,r23,r24) = (i, 1,1,1,1,1 



Thus, the task of computing the two-loop amplitude reduces to determining the coefficients 
Cj jCT . as functions of the external momenta. This is achieved by applying generalized cuts to both 



sides of eq. (4.1 ) which has the effect of turning either side into a contour integral in the complex 



plane. More specifically, the left hand side of eq. (4.1 ) is turned into a product of tree amplitudes 



whose external states coincide with either the external states of the two-loop amplitude or with 
the states traveling along the propagators that are becoming on-shell. On the right hand side of 
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eq. (4.1), the generalized cut has the effect of removing all the integrals that do not contain the 



propagators being cut, and turning the remaining integrals into contour integrals. 





19,o-i 





20,0-1 



1 21,(71 




22,01 




23,o-i 



24,oi 



Figure 7: The integral basis in terms of which the two-loop six-point amplitude is expressed. The 
integrals shown here are recorded in the o\ permutation where the external momenta are labeled 
clockwise starting with k\ at the position of the arrow. The basis is overcomplete, giving rise to 



non- uniqueness of the integral coefficients in eq. (4.1). 



To limit the number of terms remaining on the right hand of eq. (4.1) as much as possible, 
it would therefore be natural to start by imposing as many simultaneous four-dimensional cut 
constraints as possible, which at two loops would lead us to consider octacuts. Thus having 
obtained the coefficients of integrals that admit an octacut, we could then proceed to relax 
one cut constraint to allow the contributions of the double box integrals in Fig. [7] and in turn 
determine their coefficients. This, however, is not the route we will be following: instead, we will 
focus our attention on the heptacuts of eq. (4.1). Indeed, imposing seven on-shell constraints on 
two four-dimensional loop momenta leaves an unconstrained parameter z, and this makes it easy 
to survey the global poles of the integrand (to be defined properly in Section 4.1.3) on the right 
hand side of eq. (4.1). This in turn will make it straightforward to determine the tori encircling 
the global poles, which we shall refer to as leading singularity contours. 

It is important to distinguish between these leading singularity contours and the maximal- 
cut contours. Indeed, as explained in ref. |34) . the maximal-cut contours are particular linear 
combinations a of these tori whose coefficients are determined by the requirement that any 
function that integrates to zero on 'R D x ~R. D should also integrate to zero on a. This constraint 
ensures that two Feynman integrals which are equal, possibly through some non-trivial relations, 
will also have equal maximal cuts. As argued in ref. |33], multidimensional contours a satisfying 
this consistency condition are guaranteed to produce correct results for scattering amplitudes 
in any gauge theory, not only N = 4 SYM theory. It would therefore be very interesting to 
determine these maximal-cut contours. 
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Such an analysis is, however, not currently possible: the set of linear relations between two- 
loop integrals includes the set of all integration- by-parts (IBP) identities between the various 
tensor integral^ arising from the Feynman rules of gauge theory, and at present a complete 
knowledge of these relations is not available. Nevertheless, all the details needed to carry out this 
analysis (in particular, the enumeration of the global poles of the integrand and the expressions 



for the cut integrals) are provided in the intermediate results in Section 4.1 and in Appendix 
|A] and the correct heptacut contours can therefore be determined immediately once the IBP 
relations do become available. 

With the cyclic ordering of the external momenta shown in Fig. [7j the two-loop six-point heli- 
city amplitudes in J\f = 4 SYM theory admit nine different heptacuts, dictated by the eight double 
box integrals and the factorized double box integral in Fig. [7j respectively l2,a\ , ^3,01 > ^4,<xi 1 ^5,01 > 
^6,cti ? -^7,0-1 > -fi6,CTi j -^17,(ti an d Ix ai- We label these heptacuts respectively #1, . . . , #9 and study 
them in turn in Sections 14.11 and I A. 211 A. 9l 



In attempting to obtain the integral coefficients Ci :(Tj in eq. (4.1 ) from generalized four-dimensional 
cuts, one encounters two technical issues. The first point is that the basis is overcomplete, and 
the integral coefficients are therefore not uniquely defined. This feature will manifest itself as the 
appearance of free parameters in the solutions of the linear equations satisfied by the integral 
coefficients. This in turn means that one has to set a subset of the integral coefficients equal 
to specific values in order to obtain unique solutions for the remaining coefficients. The non- 
uniqueness of the integral coefficients is accounted for by the existence of various linear relations 
between the integrals in Fig. [7J as was explained carefully in ref. |73j . 

The second point is that the coefficients of the ^-integrals /i4, CTl and (thus called because 

their integrands contain factors involving the (— 2e)-dimensional part of the loop momenta) are 
of 0(e) in the dimensional regulator and hence are not obtainable from four-dimensional cuts 
|108j . As we restrict ourselves to taking four-dimensional generalized cuts in this paper, we shall 
therefore not be concerned with these integral coefficients. 

4.1 Example: evaluation of heptacut #1 

This section is intended as a pedagogical example of the use of the leading singularity method 
to obtain integral coefficients of two-loop scattering amplitudes in N = 4 SYM theory. In this 



example, we evaluate both sides of the two-loop equation (4.1) resulting after imposing the 



simplest of the nine heptacuts that the six-point amplitude admits, labeled This heptacut 



is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 



K\ = fci23 K 2 = hi K 3 = 

K 4 = h K 5 = k 6 K 6 = 0. [ ' ' 



After having evaluated this heptacut of eq. (4.1), we will set up equations satisfied by the full 



integral coefficients (i.e., including both parity-even and odd parts) and solve these equations 
explicitly. Along the way, we also present the 22 tori encircling the global poles of the integrand 
associated with this heptacut. Finally, we comment briefly on how the linear equations obtained 
from the remaining heptacuts #2, . . . , ^9 (details of which are provided in Appendix |A|) can be 
used to obtain the parity-even part of the two-loop six-point integral coefficients directly and 
report agreement with the results originally found in ref. |95| . 



4 By "tensor integral" is meant an integral whose integrand's numerator contains powers of the loop momenta 
contracted into external vectors; the corresponding standard Feynman integral with a 1 in the numerator is 
referred to as a "scalar integral" in this terminology. 
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4.1.1 Heptacut #1 of the right hand side of eq. (4.1) 



Applying heptacut #1 to the right hand side of eq. (4.1 ) will leave the linear combination of cut 
integrals shown in Fig. [8] below. 




+ C 18,CT1 



Figure 8: The integrals remaining on the right hand side of eq. (4.1 ) after applying the heptacut 
labeled $4. The cut propagators are illustrated by the inclusion of an additional orthogonal line. 
This heptacut is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 
given by eq. (4.3). 



We will use the loop momentum parametrizations in eqs. (3.12)-(3.13). Furthermore, it will 
be convenient to define the spinor ratios 

K\-k 12 -\s 12 



Pi 

Pa 



(K\k 6 ) 
(K\K\) 

\> 1 



and their parity conjugates 

2[K*k 6 ] 
2\K\K\\ 



P\ 



P2 
Ql 



1 2 



Ql 



(K"-\^ 2 \K{-) 



Pa 



2[K\K\] 



K\-k 12 -\si2 

~ (k\-\$ 12 \k\-) 

2(K\K\) ■ 



1 3 



2(K**i) 



2[^fci] 



(4.4) 



(4.5) 



This heptacut belongs to case I treated in Section [3] each of the vertical propagators in the 
double box integral l2,a l is part of some three-particle vertex. There are thus six kinematical 
solutions (shown in Fig. [2| to the on-shell constraints. The heptacut of the double box integral 
receives contributions from each of the kinematical solutions. These contributions were found in 
Section [3] (see eq. ( 3.40[ )) and take the form 



dz Ji (z) 



(4.6) 



where, using the notation in eqs. (4.3)-(4.5), the Jacobians are 



Ji(z) 



1 



327172 



x < 



(K\-\tl\K\ 
{Kl~\$\\K\ 

{k\-\$\\k\ 

{Kl\t\\K\-)z(z 



z(z 


- P')) 


-1 


for 


i 


= 2,6 


z(z 


-*>) 


-1 


for 


i 


= 4,5 


z(z 


- Qi)] 


-1 


for 


i 


= 1 


z(z 


- Qi)) 


-1 


for 


i 


= 3 . 



(4.7) 
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Accordingly, the heptacut right hand side of eq. (4.1 ), displayed in Fig. [8j receives contributions 
where the subscript in 



from each of the six kinematical solutions of the form 

c 22,cri , Ci9 jCTl 



<j> dzJi(z) ^c 2 , CTl + jj^ 



Cl8, CTl (^l +fe 6 ) 2 



(^i-fci) 2 



(4.8) 



indicates that the function is to be evaluated in the parametriza- 



tion (3.12)-(3.13) with the parameters set equal to the values 



a\ = 1 

0.2 = 



ft 
ft 



(4.9) 



and those given in Fig. [2] with the functions ft (2) and ft (2) quoted below solutions £5 and 56 
being given by 



(3 A (z) 



' 2(k 5 k 4 )(z - P 4 ) 

lk 6 k 4 ](z-Pt) 
' 2[k 5 k 4 ](z - P') ■ 



(4.10) 
(4.11) 



Instead of displaying all six contributions individually, we can make use of the fact that the 
kinematical solutions come in three parity-conjugate pairs. Namely, the contributions coming 
from two parity-conjugate solutions are obtainable from each other by the replacements 



{ij) i — ► [ij] 

{Kr\f\KJ) <— ► (K7 

Pi <— > p: 

Q l <— ► Q'i 

Q!3 i > a.4 

ft > ft 



\K7 



(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 



where the replacement rules (4.16)-(4.17) specify that if one solution has, e.g., z = 03 as the free 



parameter, then the parity-conjugate solution should be understood as having z = 04 as the free 
parameter. 

Putting everything together, the result of applying heptacut #1 to the right hand side of eq. 



(4.1) is 



1 6 / 

2^2 f dzJi(z)Ki{i 
i=l 



(4.18) 



where the Jacobians are given in eq. (4.7) and the kernels evaluated on the six kinematical 
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solutions (illustrated in Fig. [2| are 

1 C22. 



Ki{z) 
K2{z) 



C2,a 1 



C 2,tri 



C19. 



2 <1^-|# 12 |^->(P« - P«) 2 (K\-\h\Kl-)(Pl - Pg) 



(4.19) 



c 22,cti 



Cl9,eri 



2 (^-|^ 12 |^-)(z - P«) 2 (tfJ-lfclK*-)^ - P«) 
c 18 , CTl (^-|^ 6 |^-)(z - Pf) c 8 , CTl (A-M^ 6 |^-)(z - PT) 



"12 



i*r>(*-^: 



(AM^|^-)(z-P 3 -) 



a 3 (z) 

K 4 (z) 
K 5 (z) 



parity conjugate of Ai(z) (obtained by applying eqs. (4.12 1-(4.17|) 



parity conjugate of K2(z) (obtained by applying eqs. (4.12 1-(4.17|) 



(4.20) 

(4.21) 
(4.22) 



C 2,ai 



C22,. 



1 



Cl9, CT1 



2 (Kt^KtKz ~ P2) 2 (Kl-\UK\-){z - P 3 ) 
ci 8>(7l (AT^ 1 # 6 I ) (z - Pi) c 8 , ai ( X^ I # 6 I Xj" ) (z - Pi ) 



(Kt\t a \KtKz - P 2 ) (AM^|A^)(z - P 3 ) 

parity conjugate of As(z) (obtained by applying eqs. ( 4,12[ )-( |4,17 )) . 



(4.23) 
(4.24) 



4.1.2 Heptacut #1 of the left hand side of eq. (4.1) 



For MHV configurations in J\f = 4 SYM theory, the ratio of the two-loop amplitude to the 
corresponding tree-level amplitude is independent of the distribution of the helicities of the 
external states |109[ 1110) . Without loss of generality, we will therefore assume throughout the 
paper that the helicities of the external states are (l - , 2~, 3 + , 4 + , 5 + , 6 + ). In this case, the result 



of applying heptacut #1 to the left hand side of eq. (4.1) is 



W dzJ^YlAfn 

i=i 7=1 



(4.25) 



Si 



where the cut amplitude evaluated on the six different kinematical solutions yields 

' 1/1 



n-4 



tree / 
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-A 



tree 



X < 



Si 



J 3 (z) \z 

_ 1 (_L 

Ji(z) \z—, 



1 



Pi 



1 



J 5 (z) \Z-Pj 



1 



Z-Pz 

1 

Z-P3 







for i = 3 

for % = 4 

for i = 5 

for i = 1,2,6, 



(4.26) 



with the tree-level amplitude given by 

^tree 



i(12) s 



+++ 



<23)<34)(45)<56)<61) 



(4.27) 



The expressions in eq. (4.26) can be obtained by first exploiting momentum conservation to 
simplify the heptacut amplitude n| = i ^ ree | as much as possible and then substituting the 



parametrization of the loop momenta in eqs. (3.12 )-(3.13) to obtain the heptacut amplitude as 



a function of z. To further simplify, one can make use of the fact that, for N = 4 SYM theory, 
the function 

m 



<Pi(z) 



A 



tree 



Mz)l[Aj™(z) 

i=i 



(4.28) 



Si 
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only has simple poles in z and has residues ±1 or at finite poles, and at infinity |52[ I71j|j 
Denoting all finite poles of <fi(z) by Xj, these facts combined imply that 



Res z=Xj <fi{z) 



(4.29) 



It is from this latter form of ifi(z) that the expressions in eq. (4.26) were extracted. 



4.1.3 Extraction of integral coefficients 



To summarize the results of Sections 4.1.1 and 4.1.2 we find that applying heptacut #1 to both 



sides of eq. (4.1) produces the equation 



6 . 6 

J2f dzJi(z)l[Af»(z 



Si 



I f 

j^2f dzJi(z)Ki(z) 
i=i 



(4.30) 



where the Jacobians J%{z), the cut amplitude rjj=i Aj^^z) and the kernels K%{z) are given in eqs. 



(4.7), (4.26) and (4.19)-(4.24), respectively. The kernels K{{z) contain the integral coefficients, 



and by making various appropriate choices of the contours Tj, eq. (4.30) produces a system of 



linear equations which can be solved to obtain the integral coefficients. 



Before proceeding to discussing how equation (4.30) may be used to determine the integral 



coefficients of an amplitude, let us first remind ourselves of the relation of this equation to the 



original two-loop equation (4.1). As we found in Sections 4.1.1 and 4.1.2, the on-shell constraints 



in eqs. (3.3)-(3.9) and (4.3) for heptacut #1 are solved by setting the parameters on,j3j equal to 



the values quoted in eq. (4.9) and in Fig. |2| with the spinor ratios P\, Q\ etc. and the functions 
/3s(z) , ^(z) being given in eqs. (4.4)-( [4~5| ) and ( |4.10[ )- p.ll ), respectively. For any of the six 



solutions to the on-shell constraints shown in Fig. |2j one of the loop momentum parameters 
Oi,/3j is set equal to an unconstrained complex parameter z. 

The leading singularity contour is, by definition, a torus consisting of circle factors centered 
around the seven on-shell values of the parameters left fixed and where, in addition, one makes 
a choice of contour for the unconstrained degree of freedom z. The contour in z can be chosen 



to encircle the poles of the Jacobians (4.7) or the poles of the loop momentum parametrization 



(3.12)-(3.13) at which one of the loop momenta becomes infinite (see eqs. (4.10|-( 4.11| )). The 
point (ai, . . . , CK4, fa) £ C 4 x C 4 encircled by the torus is referred to as a global pole of 



the integrand of the right hand side of eq. (4.1). 



5 This property holds for MHV amplitudes (at any loop order) only; however, for N fe MHV amplitudes, the 
residues can be expressed economically in terms of superconformal i?-invariants |35| . 



23 



Defining 

r = C ai (l) x C a2 (0) x C A (0) x C A (1) , (4.31) 

the six sets of on-shell values of the loop momentum parameters quoted in Fig. [2] combined with 
the various possible choices of contours in z, give rise to the following 22 tori corresponding to 
heptacut $1, 



Ts 1 ,o = 


T 


X 


C a3 (P?) x C a4 (0) x Cp 3=z (0) x C Pi {Q) 


T Si,P? = 


T 


X 


C a3 (P'i) x C a4 (0) x Cp 3=z {PZ) x c P M 


Ts 2 ,0 = 


T 


X 


C a3= M x C a M x C%(Ql) x CpM 


t s2,p; = 


T 


X 


C a3=z (P') x C Q4 (0) x C P ,{Q\) x C A (0) , i = 1,2,3 


T s 3 ,o = 


T 


X 


C Q3 (0) x C Q4 (Pi) x Ca,(0) x C^ 4=2 (0) 


T S 3 ,Pi = 


T 


X 


C Q3 (0) x Co, (Pi) x ^(0) x C^ 4=2 (Pi) 


Ts A ,o = 


T 


X 


C Q3 (0) x C Q4=2 (0) x Cfy(0) x C A (Qi) 


T S 4 ,Pi = 


T 


X 


C a3 (0) x C a4=z (Pi) x C^(0) x C^(Qi) , % = 1,2,3 


T Ss,0 = 


T 


X 


C a3 (0) x C Q4=2 (0) x Cfyfaiz)) x (7^(0) 


T Ss,Pi = 


T 


X 


C Q3 (0) x C Q4=2 (P) x Cfa{fa(z)) xC A (0), i = 1, . . . ,4 


Ts e ,o = 


T 


X 


C Q3=2 (0) x C Q4 (0) x C^ 3 (0) x C P M{z)) 


T Se,P' = 


T 


X 


C a3=z (P') x C Q4 (0) x CpM x C^ifi^z)) , i = 1, . . . ,4 



(4.32) 



where the indices in T$ i x refer to the kinematical solution Si associated with the given torus, and 
X is the pole around which the z-contour is taken. Furthermore, C aj (X) (e.g.) denotes a small 
circle in the <x,-plane centered around atj = X; in addition, we write C aj=z {X) whenever the ay- 
variable in question is left unfixed by the heptacut constraints. Finally, the spinor ratios P\,Q\ 



etc. and the functions ^(z), (3a(z) are given in eqs. (4.3)-(4.5) and (4.10 )-(4.11 1, respectively. 

The relation of eq. ( 4,30[ ) to the original two-loop equation (4.1 ) can now easily be stated: the 
former equation is obtained from the latter by changing the integration contour from the real slice 
R D x R D to an arbitrary linear combination of the 22 tori given in eq. rt4.32|. The contributions 
from, e.g., the first two terms in this linear combination as t ,qTs 1 .o + 0'S 1 ,P'Ts 1 ,p* + • • • are found 
by integrating out all parameters except z = ^3; this leads to the Ti-integral in eq. (4.30) where 
ri = a Sl fiCf3 3=z (0) + a Sl ,p'Cf3 3=z (P'). 



Thus having explicitly stated the relation of eq. (4.30) to the original two-loop equation 



(4.1 ), let us now return to the question of how the equation (4.30) may be used to determine the 
integral coefficients of an amplitude. As explained in the beginning of this section, the class of 
contours a i,jTSi,Xj that will produce correct results for the integral coefficients in any gauge 
theory amplitude are those that annihilate all functions that integrate to zero on the real slice 
R D x R D . Determining such contours requires the knowledge of all integration- by-parts identities 
at six points; however, as a complete knowledge of all such relations is presently not available, 
we will here proceed as in ref. |73| and assume that any contour is valid. Indeed, as already 
mentioned above, the purpose of this section is mainly to provide a pedagogical exposition of the 
use of the leading singularity method to obtain integral coefficients of J\f = 4 SYM amplitudes. 

Below we will use the following notation: C e (Xk) denotes a circle centered around z = X^ 
of some appropriately small radius e (i.e., small enough to not enclose any other poles), and 
r» = 5i t jC e {Xk) denotes a contour which is zero on all six sheets except for the sheet supporting 
kinematical solution Sj\ on this sheet, the contour is a small circle centered around X^. From 



the equation (4.30) we then find, for example, that 
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setting Ti = 5i$C e (P\) produces the equation 



C-2 



C22,, 



Cl9,< 



,<7l 



2 {Kl-\$ 12 \K{-){P X -P 2 ) 2 (^-|^ 1 |iTj-)(P 1 - P 3 

-2 7172 p 1 (^-| J r 5 |^-}^ 



tree 



setting Ti = 8i$C e {P 2 ) produces the equation 



1 / 1 



C22,, 



c 18 ,a 1 (Kl-\h\K b 1 ~)(P2-Pi) 



setting Ti = Si^C^Pz) produces the equation 



1 / 1 



C19,, 



4\ 2(^-1^1^- 



c s ,a 1 (K b 2 -\h\Kl-)(P 3 -Pi) 



2 7 i72P 3 (P3 - Pi){Kl-\$ 5 \K\-)A^' 



setting Ti = 5i t2 C t (P') produces the equation 



1 



1 



C-2., 



C22 



,01 



1 



Cl9,cri 



4 V 2 ( ^-|^ 12 |^-)(p._p.) 2 { Kl\t l \K»-)(P{-P;) i 

setting Ti = 5i t2 C t (P 2 ) produces the equation 



1 



C22. 



c 18 , (T1 (^-|^|^-)(P 2 «-P 1 ') 



2(Kj-|^ 12 K 



setting Ti = 5i t2 C t (P') produces the equation 



1 / 1 



Cl9,cri 



C8)<71 (^-|^|^-)(p 3 «-p;) 



4\ 2 {K \-\^\K b 2 -) (K\-\h\K b 2 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



(4.38) 



We have here deliberately chosen an overcomplete system of equations as a consistency check on 



the method. We in fact find that the six equations (4.33)-(4.38) are consistent, and that expressed 
in terms of the quantity 



Y = -^ ll2 P,{Kl\t h \K\-)A^ ++++ (l- { S-^t\ iP l ? 



-i 



(ifn^rxp-po 



the solution of the equations (4.33)-(4.38) takes the following form 



C 2,CT1 

C22,cri 
Cl9,o-i 
Cl8,cri 

C8,cti 



Y 


2(K*r\t 1 \K*r){Pt - ps)y 



{k\-\^\k\-) y 



(4.39) 



(4.40) 
(4.41) 
(4.42) 
(4.43) 

(4.44) 
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In the following it will be useful to consider the two-loop amplitude normalized with respect to 
the tree-level amplitude, 

| (2) . _^ 

nci^haj (4.45) 



M. 



(2) 



A 



6,MHV 



6,MHV 
Atree 
MHV 



1 

4 



=1,...,24 
=1....,12 



where the normalized coefficients are defined by 



Atree 
MHV 



(4.46) 



because this object is independent of the distribution of the helicities of the external states 



(2) 



.. Moreover, the integrand of M. 



(2) 



6, MHV 



[ID21 EDI; that is, for example, M^'__ ++++ = M, 
can be decomposed into two terms that are even and odd under parity and which respectively 
coincide with its real and imaginary parts. 



Thus, the parity-even part of the coefficients found in eqs. (4.40 )-( 4.44 ) is obtained by dividing 
by the tree-level amplitude A^ff , 



and taking the real part, yielding 



Re 



Re 



Re 



Re 



Re 



c 2, 



A % 



+4 



A tree 

A +++ + 



Cl9. 



+++H 
C18,cti 



■-++++ 

c 8,cri 



/(tree 

A ++++ 



2S45S 



2S56(S123S234 ~ S^S2%) 



(4.47) 
(4.48) 
(4.49) 
(4.50) 
(4.51) 



where the equalities have been found to hold numerically. The expressions on the right hand sides 



of eqs. (4.47 )-(4.51 ) are in agreement with the results originally found in ref. [95j and reproduced 
by the leading singularity method in ref. |73j . Due to a Ward identity [109T I111] valid for M = 4 
supersymmetry, the coefficients of the cyclically permuted integrals can simply be obtained by 
cyclic permutation of the results found here. 



For the remaining heptacuts #2, . . . , #9 (details of which are provided in Appendix |A|) , the 
leading singularity method provides similar linear equations satisfied by the integral coefficients. 
However, when solving the linear equations for these more complicated heptacuts, the overcom- 
pleteness of the basis in Fig. [7] entails a problem, carefully discussed in ref. |73j : due to linear 
relations between various integrals in the basis, the integral coefficients are not unique. This 
feature will manifest itself as the appearance of free parameters in the solutions of the linear 
equations. Thus, one has to set some of the integral coefficients equal to specific values in order 
to obtain unique solutions for the remaining coefficients. 

This "gauge fixing" can be easily done for the parity-even part of the integral coefficients as 
analytic results for these were already obtained in ref. [95], and one can proceed as follows. For 
a given heptacut, encircle a pole and its parity conjugate to obtain equations of the schematic 
form 

aci + /3c 2 = 7 (4.52) 
a*a + (3*C2 = 5 (4.53) 
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where, for example, a* denotes the complex conjugate of a. From this pair of equations imme- 
diately follows 



(a - a*) cf en + ifi - /H = ilm ( ^tl J C 4 - 55 ) 

where the parity-even part of the coefficients is simply given as the real part of the (normalized) 
coefficient, 

c|™ n = Rec^.. (4.56) 



The pair of equations (4.54 )-(4.55 ) thus allows one to solve directly for the parity-even part of the 
integral coefficients. In order to obtain unique answers for these, we make the following "gauge 
choices": we set 

C U,aj — S a j (ei) s a J (12)S ( 7 j (123) 

~?V2% = s aj (456) (s CTj -(345)S CTj (456) ~ ^ (12) V; (45) ) } for j = 1, . . . , 12 (4.57) 

^1% = o 

and also recall that ci4 )CT . = Ci5 )(7 . = for the full coefficients of the //-integrals. With this 



choice, the equations produced by taking leading singularities of both sides of eq. (4.1 ) and then 



projecting out the parity-odd part of the coefficients in analogy with eqs. (4.54 )-(4.55 ) have a 
unique solution which is in agreement with the results for the parity-even coefficients originally 
found in ref. [95j; in particular, one finds that 



3g£ = • • • = = , (4.58) 

in agreement with the parity-even part of the amplitude being dual conformally invariant. 

4.2 Leading singularities vs. loop-level recursion 

In this section we report on a comparison between results for the full (i.e., parity-even and odd) 
two-loop six-point MHV J\f = 4 SYM integrand as produced by the leading singularity method 
on one hand and recent predictions in the literature based on a BCFW-like loop-level recursion 
relation |60| [63] on the other. Because the results of these papers are expressed in terms of a 
different basis from the one used in this paper, it is obviously not meaningful to check agreement 
between individual integral coefficients in the two representations. 

A quantity that can be meaningfully compared is the two-loop integrand: in general, in the 
planar limit of any field theory, the loop integrand is a well-defined rational function of the ex- 
ternal momenta (which for example can be thought of as being produced by the Feynman rules). 



In our case, the two-loop integrand is the quantity under the integral sign in eq. (4.45), obtained 
as the sum of the integrands of the 24 basis integrals Jj jCT . , weighted by the integral coefficients 
Cj jCr . and symmetry factors rj, where the summation is taken over all dihedral permutations aj 
of the external momentum labels. Agreement between the two-loop integrand as computed by 
either method would imply agreement between the integrated expressions; that is, the results for 
the amplitude]^] 

6 More accurately, the object we are comparing is the strictly four-dimensional part of the integrand. Indeed, 
as remarked above, the //-integrals , /is ><7 , give rise to contributions to the integrand which are 0(e) in the 

dimensional regulator and hence not obtainable from evaluating leading singularities in strictly four dimensions. 
Thus, our findings of agreement between the integrands should be interpreted as a statement concerning the 0(e°) 
part exclusively. 
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We have performed the comparison of the two-loop integrands numerically, by verifying agree- 
ment to high accuracy for a large number of randomly selected external and internal momenta. 
Accordingly, the remainder of this section will be devoted to discussing how, given a set of 



randomly generated momenta, one may evaluate the integrand of eq. (4.45) - that is, how the 
integral coefficients are obtained, and how the integrands of the basis integrals are added in 
a meaningful way. 

The integral coefficients are determined in analogy with the procedure explained in Section 



4.1.3 However, in that context we were concerned with obtaining analytical results for the 
coefficients and could find the coefficients of all dihedral permutations of, for example, I2 ai by 
applying the appropriate permutation to the algebraic expression for C2 ai ■ This is obviously not 
possible when one is aiming to find numerical results for the coefficients. Instead, one must act 



on the external momentum labels implicit in eq. (4.30) with each of the dihedral permutations 
(Jj £ Dq in turn so as to produce distinct linear equations for c-i^y Again, because of linear 
relations between the basis integrals, the solutions of the linear equations will, for the more 
complicated heptacuts #2, . . . , #9, contain free parameters. Accordingly, one must set some of 
the coefficients equal to specific values in order to obtain unique solutions for the remaining 



coefficients. In analogy with eq. (4.57), we choose the "gauge fixing' 



Cll.cr, = S a . ( 61 ) S a , (12) S a . (123) 

Cl2, CTj = (456) (^.(345)^ (456) - ^.(12)5^(45)) } for j = 1, . . . , 12 (4.59) 
C24, CTj = 



whereby c°f d a . = c°2 d a . = c^f^. = 0. We then find unique results for the remaining coefficients 



^odd 7;odd 7;odd 

with the property that ReCj ff) = cf^ n . 

Expressed as functions of internal and external momenta, the integrands of the basis integrals 
cannot be added in any meaningful way as the value of any term would depend on the labeling 
of the internal lines of the corresponding graph (i.e., which propagators are labeled t\ and £2)- 
To remedy this, the integrand must be expressed in terms of dual x-space coordinates, defined 
by 

1, . . . , 6 (mod 6) 
1,...,6 



Xi - 


Xi+1 


x aj (i) 


- x 7 


Xa 3 {l) 


- x 8 


3^(6) 


- x 7 


Xa,(6) 


- x 8 




X ij 




(4.60) 



7,..., 12 



i,j = 1, ... ,8 , 

with the additional requirement that, for any given graph, i\ and £2 be offset by appropriate 
translations by external momenta so that all propagators take the form -\- . Finally, the integrand 
must be symmetrized in x-j and x%. Namely, any assignment of these points to a given graph will 
fail to be invariant under vertical reflections of the graph; to ensure that the value of the integrand 
is not dependent on how its contributing graphs happen to be drawn, one must therefore average 
over the two possible assignments of these points. The integrands of the basis integrals in Fig. [7] 
have been presented in Section [A. 1| for convenience. 

In summary, given a set of random internal and external momenta, the evaluation of the 



integrand of eq. (4.45) proceeds in three steps. First, the integral coefficients are obtained by 



solving the linear equations that follow from taking the leading singularities of eq. (4.45 ). Second, 
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the integrands of the basis integrals are computed after converting the momenta into the dual x- 



space (which is achieved by solving eq. (4.60) and choosing, e.g., the base point xq = 0). Finally, 



the intermediate results are combined, weighting the contributions by the appropriate symmetry 



factors (4.2) of the integrals. This is essentially the procedure followed by our code [112J which is 
available online. As a simple consistency check of the code, we remark that the results produced 
for the integrand indeed satisfy crossing symmetry; that is, 

integrand of M 6 ^ HV (a?i, ...,x 6 ) = integrand of M^mhv( x ^(i)> • • • ' x ^(6)) for a j e D 6 ■ 

(4.61) 

We have evaluated the two-loop six-point MHV integrand for a large number of randomly selected 
rational moment^] and in all cases find agreement with |60| I63j to high numerical accuracy |96j . 
In Table 1 below we have provided a few sample points to allow the curious reader to reproduce 
our results. Further data points can be generated by the Mathematica notebook [112] available 
online. 



(xi,...,x 6 ) 



(x 7 ,x 8 ) 



integrand of M. 



(2) 



6, MHV 



2; 2> u ' u / ' V 12' 6' 4' u / ' 

4 11^ (_7 29 13 _ 1 
3' 12' 12' u / ' V 4' 36' 18' 18/ 



18' 54' 27' 



27 



), (0,0, 0,0) 



'3X1 
, 8 ' 24 ' 6 ' 6 

157 _3 _7_ 
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V) (1 
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24' 6' 
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69 _I _1Z 
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i) 



(73 _5_ _^3_ 19 \ fO 0) 
\90> 54> 135 ' 27/ ' V u ' u ' u ' u / 



1 1 n 1\ (79 _J_ 1 J5\ 
12' 3' u ' 4/ ' V96' 24' 8' 32/ 



(59 _1 J_ 
\ 48 ' 6 ' 32 ' 



24 ' 8 ' 32 / ' 

7_\ (157 _13 1 __3 
32/ ' V 96 



V224' 28' 8' 32 



24' 8' 32) ' 

), (0,0,0,0)' 



1 1 

.4 



i 2 i) 

) 4) L ) ; 

(0,0,2,0) 



(1,|,0,2), 
(0,0,1,1 



2' 3' 3' u / ' 

(i, 0,0,0 



31230 748 253 
" 22 094 130 240 



994 276 085 
981961344' 



4 777 009 838 357 
201230 662 913 280 



+ 



1802 603 853 899 ■ 
259 652 468 275 200 ' 



3 393 545 258 977 272 
16 669 297 265 



13 015 877 862 533 664 ■ 
183 362 269 915 ' 



Table 1: Values of the two- loop six-point MHV integrand (normalized with respect to the tree- 
level amplitude) at three randomly chosen sets of points (x±, . . . , xq) and (xy, xg) in dual x-space, 



respectively encoding external and internal momenta as prescribed by eq. (4.60). The parity-even 



(2) 

and odd parts of the integrand of Mg j^ H y respectively coincide with its real and imaginary parts. 
Further data points can be generated by the Mathematica notebook [112J available online. 



Finally, let us observe that we only made use of the assumption that the two-loop amplitude is 



MHV when evaluating the heptacuts of the left hand side of eq. (4.1 ). The form of the heptacuts of 



signature, the first n — 2 can be chosen 



7 To generate n rational momenta which are lightlike in (+, — , - 
as arbitrary Pythagorean quadruples (for example, generated by using the parametrization (m| + m\ + m|, m§ — 
m\ — ml, 2mim3, luiim?,) with rrn £ Z) normalized by their || • ||i-norm. To ensure that the n-th momentum will 
be lightlike and satisfy momentum conservation, the (n— l)-th momentum is obtained by generating an additional 



Pythagorean quadruple £ of unit || • |]i-norm and then rescaling it by the constant a 
k n = — (X^r^Ti 2 ^ + a i s lightlike and rational. 



' whereby 
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the right hand side of eq. (4.1 ) is independent of the external helicities, and the results presented 
in this paper can therefore straightforwardly be extended to obtain the NMHV integrand as well. 



The two-loop six-point integrand as computed in this paper was expressed in terms of the basis 
in Fig. [7] which does not include integrals containing subloops with less than four propagators. 
The exclusion of such integrals from the basis owes to the observation that cutting two propa- 
gators would factor out a subtriangle or sub-bubble — but the latter integrals are known not to 
contribute to one-loop amplitudes in J\f = 4 SYM. This strongly suggests that the uncut two-loop 
integral would appear with zero coefficient if included in the basis expansion. 

One possibility which cannot be rigorously ruled out by this argument, however, occurs when 
the two-particle cut is shared between several integrals containing subtriangles or sub-bubbles: 
in principle, the coefficients of the respective cut two-loop integrals could be nonzero, but such 
that the contributions cancel. 

Ruling out such a scenario completely would require extending the analysis of this paper to 
consider hexacuts, pentacuts etc. Two-loop integrals whose subloops contain at least three prop- 
agators all admit T 8 -integration contours (analogous to those in eq. (4.32)) that appropriately 
define such cuts. Moreover, the cut integrand is a holomorphic function, and the contour inte- 
grations can thus be performed directly by means of the global residue theorem |52j . In contrast, 
two-loop integrals with bubble-subloops do not admit T 8 -contours: for example, the bubble-box 
integral rather admits a T 6 x S 2 contour. Several approaches are available to deal with bubble- 
inherited S 2 contours, among the more elegant ones is that of Mastrolia |1 13] . exploiting Stokes' 
theorem; a related approach is that of Arkani-Hamed et al. in ref. |114| . However, we leave such 
extensions for future work. 



5 Conclusions 

In this paper we have provided a check that recent results in the literature |60[ [63] for the 
full (i.e., parity-even and odd) two-loop six-point MHV integrand of M = 4 SYM theory can be 
reproduced by the leading singularity method. Equivalently, assuming the validity of refs. |60U63] . 
one can view the analysis carried out in this paper as a check that the leading singularities of 
the A/ - = 4 SYM integrand evaluated in strictly four dimensions (as opposed to in D = 4 — 2e 
dimensions) are sufficient to detect the parity-odd part. This has already been shown to be the 
case for the two- loop five-gluon J\f = 4 SYM amplitude in ref. [72j, but the six-gluon MHV 
amplitude provides a much richer testing ground owing to the larger number of global poles of 
the integrand and of integrals in terms of which the amplitude is expressed. As the main part 
of the results presented in this paper are independent of the helicities of the six external states, 
the NMHV integrand can be straightforwardly obtained by supplementing the requisite helicity 
dependent data. We leave this as an open problem. 

For any two-loop integral, the leading singularities are obtained by changing the integration 
range from IR^ x ~K D into tori of real dimension 8 (embedded in C 4 x C 4 ) that encircle the global 
poles of the integrand. As explained in ref. |34j . the maximal cuts are particular linear combi- 
nations a of these tori whose coefficients are determined by the requirement that any function 
that integrates to zero on M. D x M. D should also integrate to zero on a. This constraint ensures 
that two Feynman integrals which are equal, possibly through some non-trivial relations, will 
also have identical maximal cuts. As argued in ref. |34j . multidimensional contours a satisfying 
this consistency condition are guaranteed to produce correct results for scattering amplitudes in 
any gauge theory, not only J\f = 4 SYM theory. 
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The set of linear relations between two-loop integrals includes the set of all integration-by- 
parts (IBP) identities between the various tensor integrals arising from the Feynman rules of 
gauge theory; however, at present a complete knowledge of such relations is not available. Thus, 
an interesting open problem to be pursued once all necessary IBP relations do become available 
is to determine the maximal-cut contours that allow the extraction of integral coefficients in any 
two-loop six-point gauge theory amplitude. We expect the intermediate results provided in this 
paper to greatly facilitate this task. 
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A Details of heptacuts #2, . . . , #9 



A.l Integrands of basis integrals in dual coordinates 

The expressions below are the four-dimensional integrands of the basis integrals in Fig. [7j ex- 
pressed in dual x-space coordinates (related to the internal and external momenta through eq. 
(4.60)). The results are recorded in the o\ permutation of the external momentum labels; the 
integrands for the remaining dihedral permutations may be obtained by applying eq. (A. 2 ) below. 



integrand 1(Tl 
integrand 2 ai 
integrand 3(T1 
integrand 4(Tl 
integrand 5 (Tl 
integrand 6(Tl 
integrand? (Tl 
integrandg (T1 
integrand 9(T1 
integrand 10(Tl 
integrand 11 
integrand 12(T1 
mtegrand 13(T1 

integrand 14(T1 
integrand 15(Tl 
integrand 16(T1 

integrand 17 (T1 

integrand 18(T1 

integrand 19 CTl 

integrand 20(Tl 

integrand 2lCTl 
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The integrands for the remaining dihedral permutations <jj G I?6 can be obtained from eq. ( A.l ) 
by applying 



integrandj j(Tj {x\ , x 2 , 



,x 6 ;x 7 ,x 8 J 



integrandj ai 
integrand 



(^■(1)) X cr j (2) ) 



, x <j,(6); x 7, x 8) 
for j = l,. 

,X aj ^)]X 7 ,X 8 ) 

for j = 7,. 



12 . 

(A.2) 
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A.2 Heptacut #2 



This heptacut is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 

(A.3) 



K 2 = h 
K 5 = k 6 



K 3 = 
K 6 = 0. 



Applying this heptacut to the right hand side of eq. (4.1 ) leaves the following linear combination 
of cut integrals 



C3,c 



1 — 1 — 1 


i — i— - 


tr 
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C 2 0,a 4 fi 
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(/l + fc 6 ) 2 
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*1 *6 
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1+f )^-fc 56 -|8 56 
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P 3 
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(l+f )[^fcl] 
2[i^fci] 



(A.4) 



33 



and their parity conjugates 



[Kike] 



1 2[K\k 6 



Ql 



2" 

\K\K\] 
2\K\K\ 

S4 



1+ 



12 



{K\K\) 



2(K\K\) 



P* 

r 2 



Ql 



2\K\k x \ ' 



(l+fj^fc56-^56 



P* 



Ql 



' (k\-\^\k\-) 



2(K\k 1 ) 



(A.5) 



This heptacut belongs to case I treated in Section 3.1 and there are thus six kinematical solutions 
(shown in Fig. [2j. Parametrizing the loop momenta according to eqs. (3.12 )-(3.13 ), the on-shell 



constraints in eqs. (3.3)-(3.9) are solved by setting the parameters equal to the values 





OL\ = 1 , 
a 2 = , 



Pi 

02 



1 + 



^1 

72 



(A.6) 



and those given in Fig. [2] with 



1 + 



SA {k\k\)(z - Pi) 



(A.7) 



72 J 2{KlKl){z-P A ) 

for kinematical solution £5. The heptacut double box integral I^^i is Si=i §v .dzJi(z) where 



Hz) 



1 



327172 



X < 



1 + 1) (K\-\t\\K' 2 



Z[Z 



1 + 1) ^K\-\t\\K\-) z{z - Pi 



{K\-\t\\K\-) z(z - Of] 

(jfS-|jt*|its->z(z-oi: 



for z = 2,6 

for i = 4, 5 

for i = 1 

for z = 3 . 



(A.8) 



A. 2.1 Heptacut #2 of the right hand side of eq. (4.1) 



The result of applying heptacut #2 to the right hand side of eq. (4.1) is 



1 6 / 

7^2 f dzJi(z)Ki{z) 
i=i 



(A.9) 



where the kernels evaluated on the six kinematical solutions are 
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K 2 (z) = c 3 , CTl + 



1 
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where Ps(z) is given in eq. (A. 7) 



A. 2. 2 Heptacut #2 of the left hand side of eq. (4.1) 



The result of applying heptacut ^2 to the left hand side of eq. (4.1) is 



:J2<f dzJi(z)f[A$°°( 

i=l JTi 7=1 



(A.16) 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the six different kinematical solutions yields 

- A!!! 6 ^^ x I _^J\ for i = 4 (A. 17) 

for z = 1,2,5,6. 
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A.3 Heptacut #3 



This heptacut is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 

(A.18) 
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Applying this heptacut to the right hand side of eq. (4.1 ) leaves the following linear combination 
of cut integrals 
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and their parity conjugates 
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This heptacut belongs to case I treated in Section 3.1 and there are thus six kinematical solutions 



(shown in Fig. |2j. Parametrizing the loop momenta according to eqs. (3.12 )-(3.13 ), the on-shell 

(A.21) 



constraints in eqs. (|3.3|)-(|3.9|) are solved by setting the parameters equal to the values 
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for kinematical solution S§ . The heptacut double box integral is Ylt=i §v .dzJi(z) where 
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A. 3.1 Heptacut #3 of the right hand side of eq. (4.1) 



The result of applying heptacut ^3 to the right hand side of eq. (4.1) is 
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where the kernels evaluated on the six kinematical solutions are 
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(A.29) 
(A.30) 



where f3%{z) is given in eq. (A.22). 



A.3.2 Heptacut #3 of the left hand side of eq. (4.1) 



The result of applying heptacut ^3 to the left hand side of eq. (4.1) is 



6 . 6 
J2f dzJi(z)l[Af 

i=l 1=1 



(A.31) 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the six different kinematical solutions yields 



n-4 



tree / 



1 



16 



A 



tree 



1 



J 5 {z) \z-P 2 




z-P A 



for i = 5 

for i = 1,2, 3, 4, 6. 



(A.32) 
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A.4 Heptacut #4 



This heptacut is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 

(A.33) 



Ki = fci 
K 4 = k 34 



K 2 = k 2 
K 5 = k 5 



K 3 = 
K 6 = k 6 . 



Applying this heptacut to the right hand side of eq. (4.1 ) leaves the following linear combination 
of cut integrals 





1 


i 




h 



5 f 6 k. 



+ C13 
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We define the spinor ratios 

Pi 
Pi 

Qi 
Q 4 



l+^)Kl-(Kl+k 6 )+Kl-k 6 

<^-|(l+f )t\+h\K\-) 
(K\-\t\+^\Kl-) 
2{K\k\)[K\k\\ ' 

l+^)Kl-{Kl+k 6 )+K\-k 6 
{Kl-\$\+f & \K\-) 



p 2 
Q2 



(K>-\h\Kl~) ' 

^■^456+1*456 

(k" 2 -\^ 6 \k\-) 



l+f)^-fc4 5 -|s 45 

(Ki-\h 5 \K b 4 ~) 



ft = " 



K^kse+^sse 



{K\-\t 6 \K\-) 



(K b 5 -\^\Kl-) ' 
and their parity conjugates 



(A.34) 



P'l 



P* 

1 1 



Ql 
Ql 



1+ 



A - 4 



)^-(x5+fc 6 )+i^f-fc 6 



(A*-ir 1+ jj 6 i*rn 



(Jf*-[*-i+JJ 6 |Jf*-> 



P* 

J 2 



P* 



Ql 



K\-k & 

' (k\-\^\k\-) ' 

^■^456+ 5^456 
(l+f|)^.fc4 5 -|s 45 

<jti-ijt«i^~> 



^■fc 5 6 + gS56 
<^-|^ 56 |^-> 



(^-|^ 6 |^-> 



(K-ri^i^r> 



(A.35) 

This heptacut belongs to case I treated in Section 3.1 and there are thus six kinematical solutions 
(shown in Fig. [2j. Parametrizing the loop momenta according to eqs. (3.12 )-(3.13 ), the on-shell 
constraints in eqs. (3.3)-(3.9) are solved by setting the parameters equal to the values 





"1 = 1, 
a-2 = , 



1 + 



Si 

12 



and those given in Fig. [2] with 



Ql{P2-P±){z-P x ) 



(A.36) 



(A.37) 



(P 2 -P x ){z-P A ) 

for kinematical solution £5. The heptacut double box integral /5 j(71 is Y^a=i §t -dzJi(z) where 



Hz) 



1 



327172 



x < 



K\- 



1 + 
1 + 



$5 + ^6 



+ , 



{K\-\t\ + h\Kl 



ziz 



K\- 

Q\) 



{Kl-\$\ + UK\-)z{z-Q{) 



ziz 



p-: 



z(z - ft) 

-1 

-1 



for i = 2, 6 

for i = 4, 5 

for i = l 

for i = 3 . 



(A.38) 



A. 4.1 Heptacut #4 of the right hand side of eq. (4.1) 



The result of applying heptacut ^4 to the right hand side of eq. (4.1) is 



1 6 / 

7^2 f dzJi(z)Ki{z) 
i=i 



(A.39) 
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where the kernels evaluated on the six kinematical solutions are 



K\{z) = c 5jCTl 



1 



c 23,ct 3 



c 21,<r 3 



C20.. 



C22,, 



2(^-|^|^->(P--P*) 2^-|^|^-)( ar _Q.) 

1 C 2 4^(^-|fe 6 |^-)(P 1 '-P 3 ') 

2 <^-|^|^-)(^-|^ 45 |4-)(P 1 * - P«)(z - Q5) 



l C2 4, g6 (^-|fe 1 |^)(z-QI) + c 24 , CT8 (^-|^ 456 |^-)(P 1 »-P^) + C24 , gll <^-|^ 6 |^-)(z-Qg 
2 (^-|^|^-)(^-|^ 45 |A^-)(P 1 - - P')(z - Q«) 



(Kt\h\Kt)(Kl\h 5 \Kt)(Pi PSKz 05) 

cl3, CT e(^^l^456|g^)(g^lfell^^)(A , - gncg - q» 

(XM^ 6 |X^)(XM^ 5 |^-}(P 1 - - P*)(z - Q2) 

Cl2,. 3 (AM^ 456 |^-)(gMfe|^)(P 1 *-P5)(^-Q3) 

(^-|^ 6 |^-)<ifM^ 45 |^-)(P 1 - - P 2 *)(z - Q5) 



c 9 , CT6 (^-|fe 6 |^-)(P 1 »-P 3 ') 
(K{-\h\Kl){Pt - PZ) 

c 18 ^ 3 (Kl\hi\Kt)(z-Q'i) 
{Kt\UK\-){z-QD 



c 12 , CT8 (^-|^ 56 |^->(^-|^ 1 |^-)(P 1 * - P*)(z - Q«) cn^^-l^l^-)^ - Ql) 



(Kt\h\Kl)(K-\h 5 \Kl-)(P; - P' 2 ){z Ql) 



{KUh,\Kl-){z-QD 



(A.40) 



Ki{z) = c 5)l 



C23,, 



C21,, 



4(A-J-|^|^-)(^-|^ 45 |^-)(z-P*)(Q*-Q*) 2(^-|^ 45 |^-)(Q._Q.) 



+ 



C20,, 



c 22 ,, 



C9 ^(AM^ 56 |^-)(z-P 3 *) 



2(]fMfeK)(z-?') 2(^-|^|^->(Q--Q 3 ) (JT{-|^|^-)(z-P 2 «) 

1 e 24 , CT3 (^-|^ 6 |^-)(z - P») + c 2i , a8 (Kl$ i56 \K b 2 -)(z - P$) + c 2i , ail (I<l-\^\K b 5 -)(Ql - g|) 

2 (AM^|^-)(^-|^ 45 |^-)(z - P-)(Q- - Q5) 

ci3, CT3 (^-|^ 56 |^-)(^-|^|A^)(z - P 3 ')(QI - Qg) c 11 , g3 (^-|^ 6 |^-)(Q; - Ql) 
(K\-\UK\-)(K\-\^\Kr)(z - PZ){Q\ - Ql) (Kl\h 5 \K b 5 -)(Ql-Q' 2 ) 

c l2 ^{K\-\^\K\^{Kl%\Kl-)(z - P 5 «)(Q? - g|) 
(A 1 b -|^ 6 |^-)(^-|^ 45 |A^)(z - P 2 -)(Q- - Q5) 



(A.41) 



Kq(z) — parity conjugate of K\{z) (obtained by applying eqs. (4.12 l-(4.17l) 



K±{z) — parity conjugate of K 2 (z) (obtained by applying eqs. (4.12 1-( 4.17| ) 



(A.42) 
(A.43) 
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K 5 (z) = c 5jCTl 
1 



C 20,(T 6 



2 (^-|^ 45 |^-)(/? 3 (z) - Q2) 2 (^-i&ijtJ-x* - P 2 ) 



C22,, 



C 2 3,, 



2 (xM^ 6 |^-)(/3 3 (z) - Q5) 4 (^-|^ 6 |^-)(^-|^ 45 |^-)(z - P 2 )(/3 3 (z) - Q5) 
Cl0 ^(^-|fe|^-)(/? 3 W-QI) , c 9 , CT6 (^-|^ 56 |^-)(z-P 3 ) 



+ 



Cl8 , g3 (^-|^ 61 |^-)(/33(^)-Q?) 1 
(JCS-|jt 6 |J^-)(/3 3 (z)-QS) 



«|| 6 |^-)(z-P 2 ) 



C24 , CT 3(^-|fe 6 |^-)(^-P 3 ) 



(^-|^ 6 |^-)(^-|^ 45 |^-)(z - P 2 )(/3 3 (z) - Q5) 



C24,. 6 (^-|fe|^-)(/3 3 W-Q!) 



c 11 , g3 (^-|^ 6 |^-)(/3 3 (z)-Qg) 

(^-|^|^-)(^M^ 45 |^-)(z - P 2 )(/? 3 (z) - QJ) (X*- |^ 45 |^-)(/3 3 (z) - Q5) 

C24, CT8 (^-|^ 456 |^-)(z-P 5 ) 



2 {Kl-\UK{-){Kl-\^\Kt){z P 2 )(Mz) Q') 

1 c 2 4, ffll (^-|^|^-)(^ 3 (-z)-Q5) 

2 (^-|^|^-)(xM^ 45 |^-)(z - P 2 )(/3 3 (z) - Q5) 
c 13 , CT3 ( 1 # 56 1 Xj- ) (^M ^ 6 1 ) (z - P 3 ) (/3 3 (z) - Q* ) 



(xM^ 6 |^-)(xM^ 45 |^-)(z 


-P 2 )(/3 3 (z)-^ 


>S) 


c 13 , ff8 ^-|^ 56 |^-)(^-|^ 61 |^ 


-)(z-P 5 )(/3 3 (z 


-Q'i) 


(xM^ 6 l^->(xM^ 45 l^-)(, 


- - P 2 )((3 3 (z) - 


35) 


Cl2,«r,<^-|*«6l^l"><«4"lfel«6" 


)(z-P 5 )(/3 3 (z) 


-Q5) 


(xM^I^->(^-|^ 45 l^->(z 


-P 2 )(/3 3 (z)-C 


25) 


Cl 2, CT8 (xM^ 6 |^-}(xM^ 61 |^- 


)(z-P 3 )(/3 3 (z) 


-Ql) 


(^-|^|Xj-)(^-|^ 45 |^-)(z 


-P 2 )(/3 3 (z)-C 


25) 



X 6 (z) = parity conjugate of ^(z) (obtained by applying eqs. (4.12 1-( 4.17 1) 



(A.44) 
(A.45) 



where j3^(z) is given in eq. (A.37). 



A.4.2 Heptacut #4 of the left hand side of eq. (4.1) 



The result of applying heptacut #4 to the left hand side of eq. (4.1) is 



6 . 6 

:J2f dzMz)l[Af»( 

i=l 1=1 



(A.46) 



Si 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the six different kinematical solutions yields 



i=l 



^itree 

16 



x < 



l 

Ji(z) \z Z — P 2 



1 / 1 



1 



Je(z) \z 




for i = 4 
for z = 6 
for i = 1, 2, 3, 5 . 



(A.47) 
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A. 5 Heptacut #5 



This heptacut is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 

(A.48) 



K\ = fcl2 
K4 = /u4 



K 2 = k 3 
K 5 = k 5 



K 3 = 
K 6 = k 6 . 



Applying this heptacut to the right hand side of eq. (4.1 ) leaves the following linear combination 
of cut integrals 



c 6,< 



1 — 1 — 


1 — 1 — 


tr 

1 


— 1— 

h 



k 



C19 




Cl3, CT1 



k, k 6 

(h + k 56 f(l 2 + fc 612 ) 



(h + k 6 ) 2 (l 2 + fc 61 ) 2 



k, k 6 

(h + k m ) 2 (l 2 + fc 61 ) 2 




Cll 



01 




C 2 4,, 




+ 024,0-7 




+ C24,ctio : 



{h + he) 




[h + kei) 
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We define the spinor ratios 
Pa 



Kl-kQ + {K\+k 6 )-K\ 
{K b 2 -\Hc\+h\K\-) 



Qi 



2(K>k6) ' 
K\-k 6 +Kl-{K\+k & ) 



(K b 5 -\t\+h\K b 4 -) 
and their parity conjugates 



P* 



Ql 



Kl-ks+{Kl+k 6 )-K\ 

' (K\-\t\+h\K\-) 
[Kjk 5 ] 
2\K\k h \ ' 



' (K\-\t\+^\Kl-) 



P-2 
P5 

Q2 



P* 



Ql 



(k\-\$\+^\kI-) 

2(K\K\)[KlK\} 



2[K\k fS ] ' 



[gjfci] 
2[A**i] 



ri+jt 6 |K*-> 



2[i^](K^> 



Pa 



Qs 



p* 



2{K\k 6 ) 



(A.49) 



^5-fc61 + §S61 
'(^-[^ 61 [^-) 



~2[i^fc 6 ] 



<^-l* 6 llJ^-> 



(A.50) 



This heptacut belongs to case I treated in Section 3.1 and there are thus six kinematical solutions 
(shown in Fig. [2j. Parametrizing the loop momenta according to eqs. (3.12)-(3.13), the on-shell 
constraints in eqs. (3.3)-(3.9) are solved by setting the parameters equal to the values 



ct\ = 1 
«2 = 



01 

02 



and those given in Fig. [2] with 



Q«(P 3 - P 5 )(z - Pi] 



(A.51) 



(A.52) 



(P 3 - P x )(z - P 5 

for kinematical solution 1S5. The heptacut double box integral Iq^i is Si=i §v .dzJi(z) where 

{K»-\F 5 + t 6 \K*-)z(z-p; 

'{K\-\t\ + UK\-)*{z-Pi. 
'{Kl\t\ + UKt)z(z-Ql) 

\Kl-\t\ + UKt)z(z-Qi) 



Ji(z) 



1 
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x < 



for i = 2, 6 

for i = 4, 5 

for i = 1 

for j = 3 . 



(A.53) 



A. 5.1 Heptacut #5 of the right hand side of eq. (4.1) 



The result of applying heptacut #5 to the right hand side of eq. (4.1) is 



1 6 / 

7^2 f dzJi(z)Ki{z) 
i=i 



(A.54) 



where the kernels evaluated on the six kinematical solutions are 
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Ki(z) = C6, ffl - 



1 



C21,<7i 



Cl9 : cr s 
b — i a£ i tA-\ 



2 (^-|^|^->(P- - P-) 2 (^-|^ 6 |^-)(p. - P- 



C20,, 



C23,, 



2 - Q-) 4 (^-|^|^-)(^-|^|/^-)(p- - P'){ Z Q-) 

1 c 24 , gl (^-|fe 12 |^-)(z - gi) + c 24 ^ 4 (^ 1 b -|^ 6 |^-)(P 1 » - P 3 «) + c 24 , gl0 (^-|^ 61 |^)(z - Qg 

2 {K\-\UKl-){Kl-\UKt){Pl-P^){z-Ql) 

c 13 ^Kl-\h\Kt)(Kl~\hi\K b 5 -)(P? - P;)(z - Q' 3 ) , c 9 , CT4 (^-|fe 2 |^-)(z-Q!) 



(if 1 b -|^ 1 |^-)(^-|^|^-)(P 1 --P 2 -)(z-Q^) 



(Kt\hm~)(*-Q m 2) 



c^AKl\h\K b 2 -)(K b 4 -\$ 612 \Kl-)(p; - p 8 «)(*-q;) Cll , CT10 (x 1 b -|^ 6 |^-)(p 1 * - pj) 



(Kl-\^\Kt)(K-\h\Kt)(P' - P!)(* - 05) 



(AM^-XP'-P-) 



(A.55) 



^ / \ _ c _ 1 c 21,ctio , 1 c 19,ct 6 , 1 C 20,CT 4 

6 ' CT1 2(^-|^|^-)( z _p.) 2 {K [-\lt 6 \Kl)(z -PJ) Z{Ki-\% 6 \Ki-)(Ql-Q- 2 ) 

4 (^-|^|^-)^-|^ 6 |^-)( 2 - P-)(Q- - Q-) (^-|^ 6 |^-)(z - PS) 

1 c 24 , CT4 (A-r|^ 6 |^-)(z- P») + c 2i ^(Kl-\^ 6 \K b 2 -)(z - gT) +C24, gl0 (^-|^ 61 |^-)(Q; - g| 

2 (^-|^|^-)(^-|^ 6 |^-)( 2 - P 2 ')(QJ - Q5) 

ci3^(A 1 b -|^|^-)(^-|^ 61 |^-)(z - P 3 ')(QI - Qg) Cl0 , glo (^-|^ 56 |^-)(z - Pf) 
(Kt\h\K b 2 -)(Kt\^\Kl-)(z P a -)(Q; - Q5) (AM^I^Xz - P«) 

ci2^(A 1 b -|^ 56 |Ar)(^-|^ 1 |^-)(z - PfJCQi - Q5) Cll ^ 10 (^-|^ 6 |^-)(z - P 3 ') 



56 1 

{Kt\h\Kt){Kt\h\KfK* P')(Q'i Q5) 



(A.56) 



A 3 (z) = parity conjugate of Ai(z) (obtained by applying eqs. ( 4. 12 )-( 4. 17 1 ) 



(A.57) 



Af 4 (z) = parity conjugate of K 2 (z) (obtained by applying eqs. (4.12 )-( 4.17 1 ) 



(A.58) 
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K 5 (z) = C6, ffl - 

1 



1 



c 21,crio 



1 



Cl9,, 



1 



C20,ct 4 



2 - P 2 ) 2 - P 3 ) 2 (Kt\$ 6 \Kl)(0 3 (z) Q5) 

c 23 ^ , c 8 , CT6 <^-|^ 56 |^-}( 2 r-P 1 ) 



4 (^-|^|^->(xM^|^-)(z - P 2 )Q3 3 (z) - Q-) (^-|^ 6 |^-)(z - P 3 ) 

C 9 ^(Kl-\hl2\Kt)(Mz)-Ql) 1 C 2 4, CTl (^4"lfel2l^5~)(/33^)-QT) 

(ifM^|^-)(/3 3 (z) - QVl 2 - P 2 )(/3 3 (z) - Q«) 

l c 24 , ff4 (^-|^ 6 |/^-)(z-P 3 ) + c 24 , gT (^-|^ 56 |^-)(z-Pi) + c 24 , CTl0 (^-|^ 1 |^-)(/3 3 (z)-Q^) 
2 (^-|^K->(^-|^|^-)(z-P 2 )(/3 3 (z)-Q-) 

- P 2 )(/3 3 (z) - Q5) - P 2 ) 

c 13 ,^ 2 -\h\Kl-){Ki-\hi\Kl)(z - P a )(J3 3 {z) - Q5) Cll ^ a {K\-\h\K[-){z - P 3 ) 



(ifM^|^-)(^M^ 6 |^-)(z - P 2 )(/3 3 (z) - Q5) 

Cl2, gl fej^feK^l^l l*5~) (* - gOOW - g|) 

(Kt\h\Kl)(Kl\ft 6 \Kl)(z~P 2 )(l3 3 (z)-Q-) 

c 12 , CT7 (^-|^|^-)(^-|^ 612 |^-)(z-P 3 )(/3 3 (z)-QI) 
|^ |^-) % \K\-) (z - P 2 )(/3 3 (z) - Q5) 



(KM^|^-)(z-P 2 ) 



(A.59) 



i^^) = parity conjugate of ^5(2) (obtained by applying eqs. (4.12 1-( 4.17| ) 



where /3s(z) is given in eq. (A.52). 



A. 5. 2 Heptacut #5 of the left hand side of eq. (4.1) 



The result of applying heptacut ^5 to the left hand side of eq. (4.1) is 



W dzMz)l[Af~( 

i=i jTi 7=1 



(A.60) 



(A.61) 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the six different kinematical solutions yields 



tree / 



1 tree 



16 



x < 



1 / 1 



1 



J 4 (» \z-P 2 z-P A 





for i = 4 

for i = 1, 2, 3, 5, 6 . 



(A.62) 
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A.6 Heptacut #6 



This heptacut is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 

(A.63) 



K x = k x 

K4 = &4 



K 2 = k 2 
K 5 = k 5 



K 3 = k 3 
Kfi = fa 



Applying this heptacut to the right hand side of eq. (4.1 ) leaves the following linear combination 
of cut integrals 




{h + hvi) 



c 24,ct 3 




Cl3 





(h + k >6 ) 2 



(h + k 56 ) 2 (l 2 + k 6 ) 2 



Cl3. 



C] 




Cl3 





(h+k i56 ) 2 (l 2 + k 61 ) 2 



k, k 6 

(h+h 6 ) 2 (l 2 + k m ) 2 
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+ Ci2,, 




k 3 k 4 k 3 k 2 k 5 k 6 

(h + k 6 f(l 2 + k ei2 f (h + h, 6 ) 2 (l 2 + k 6 ) 2 (h + h 6 ) 2 (l 2 + k 61 f 



+ Cll,o 




(h + hf 




+ Cu, 



(k + k 6 ) 2 




C24 




(*i + fee) 




C 2 4,, 



We define the spinor ratios 
Pi 



Pa 

Qi 
Qa 



Kl-{K\+k 6 )+K\-k 6 

' {K\-\t\+h\K\~) ' 
(K\H) 
2(K^k 5 ) ' 



cu,, 




('2 + fc 6 12) 2 



C24 



0-7 



C24,, 




+ c 24,ctio 




('2 + fcei) 



C24,, 



^■^123-5^123 




{K>-\t 123 \K»-) ' 
p (^-|jt b 1+ )6 6 |^-) 
5 2[i^](Kpf|) ' 



^■(^+fc 6 )+^-fc 6 



<jr*-|tfi+j6 6 |X*-> ' 



Q 2 
Q 5 



-^5-^345- 1 S345 
{K\-\t\+hK~) 

2\K\K\]{KlKl) 



P 3 
P6 

Q 3 



(h + k (i ) 



' 2(K\kf.) 

(K\-\t\+h\Kt) 
2(K\k\)\K\k\\ 

\K\K\ 
2[K\k 6 ] 



(A.64) 
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and their parity conjugates 



Pi 



P* 



Ql 
Ql 



Kl-(K\+k 6 )+K\-k 6 



2[K»k 5 ] 



" (^-|it b 1+ ^|^-> 
2(K\k 1 ) ' 



*2 



1 5 



'(ArJ-|^ 123 |^-> ' 

(K*-|X*i+jt 6 [J^-) 
2(Jt*Ki>[XjJSS] 

-^"5-^345 — |«345 
^-|^ 345 |^-> ' 



P* 



Ql 



{Kjk 6 \ 
2[Kpe] 

(Jg*-[jti+jt„|Jg*-) 
2(S) 



(A.65) 



2(K\K\)[KlKl] 



This heptacut belongs to case II treated in Section 3.1, and there are thus four kinematical 
solutions (shown in Fig. |3j). Parametrizing the loop momenta according to eqs. (3.12 )-(3.13 ), the 



on-shell constraints in eqs. (3.3)-(3.9) are solved by setting the parameters equal to the values 

(A.66) 



Oil 
OL2 



1 





01 

02 



and those given in Fig. j3j The heptacut double box integral is Yli=i §r ^ z ^i( z ) wnere 

(K^\^ 1 + ^\Kl)z(z-Ql) 

[{Kl-\t\ + h\Kl-)^-Qi) / 
[{Kl^ + h^n^-Pi)^ 
\K\-\Hc\ + UKt)z{z-Pl) 



Ji(z) 



1 

327172 



x < 



for i 

for i 

for % 

for i 



1 

2 
3 

4. 



(A.67) 



A. 6.1 Heptacut #6 of the right hand side of eq. (4.1) 



The result of applying heptacut #6 to the right hand side of eq. (4.1) is 

dzJi(z)Ki(z) 

where the kernels evaluated on the four kinematical solutions are 



l -± 

4^ jr 

i=l JL * 



(A.68) 



Ki{z) 



C7,, 



C23,- 



C21,. 



4 (Kl\!t 123 \Kl-)(Kl-\fr 6 \Kl-)(a 3 (z) P})(z - Ql) 2 (^-|# 123 |i^)(a 3 (z) - PJ) 

1 C 23 : <T 3 1 C 21,cr 3 

4(^-|^ 6 |^-)(^-|^ 45 |^-)(a 3 (z)-P*)(z-Q') " 2(^-|^ 345 |^-)(z-Q5) 



+ 



1 



C21.C78 



1 



C 21,<r 7 



Cl o )fflo <^-|^ 6 i^->(a 3 ^)-pr) 



2 (ifJ-|^|X*-)(a 8 (*) - PJ) 2 (*T |*el*6 >(* - 05) 1^23^2 > - PS) 

1 c 24 , gl (^4~ 1^612 l^s") - gj) + c 2 4, g4 (^-|^ 6 |^-)(a 3 (z) - P?) 

2 {K\-\^\K\-){Kl\UKl-){a a {z) - f?)(z - Ql) 

1 c 2 ^ 7 (Kt\h e \Kt)(Mz) ~ p i) + C24, gl0 (^~|^ 6 il^~)(^ - gj) 

2 (^-|^ 123 |^-)(^M^|^-)(a 3 (z) - P«)(z - Q«) 

do^^-l^ll^X^-QI) l C 2 4 ^(^lfe6l^2"X«3W-Pl , )+C24,. 6 (^|fe|^)(^-Ql) 



(K b 4 -\t 3i5 \K b 5 -)(z-Q') 



(Kl\^ 6 \K b 2 -)(K b 4 -\^ 45 \Kl-)(a 3 (z)~P 3 ')(z 



1 c 24 , g8 (^-|^ 56 |^-)(a 3 (z) - P 2 *) + C24, gll (gM^|^-)(^ ~ gf) 

2 <|WO«Hf 345 |^-)(a 3 (z) - P«)(z - Ql) 
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c 10><T6 {Kl\# 5(i \K»-)(a 3 (z)-P;) c 13 ^Kl-\^ m \Kl}(Kl\^ 12 \Kl-)(a 3 (z)-P 1 ')(z-Q' 2 ) 



(K\-\UK-){a 3 {z)-PZ) 



{K\-\h 2 ,\K b 2 -){Kl-\UKl)(a,{z) P')(z Q') 

4- 



c w<a7 {Kt\% &1 \Kl-){z - Ql) c 13 ^Kl-\h 6 \K b 2 -)(Kl\^\Kl-)(a 3 (z) - P?){z - Ql) 



(^-|^ 6 |^-)(z-Q-) 



(xM4|^-)(^M^ 45 |^)(a3(z)-P 3 , )(^-Q2) 



Cl3 , tT10 (^-|^|^-)(^-|^ 1 |^-)(a3(^) - g?Xg - QV) 

{K\-\^\K\-){Kl\UKl-){a 3 {z) - P 2 ')(* - Q«) 

Cl3, g8 (^~l^ 5 6l^)(^4~l^6ll^5~)(«3^)--P?)(^-Q?) 

(jcM^I^-XJrM^I^-XasW - - Q5) 

cii,a 10 (Al- \h \K\~) (a 3 (z) - P 3 ') c 12iCTl (^-| jt B8 1 JH 61 \K\-)(a 3 (z) - P?){z - Q\) 



(Kt\h23\Kl)(a 3 (z) - P 2 ') 



Cll , a3 {Ki-\^\KtKz-Q' 3 ) Cl2 ^ 10 (^-|^ 6 |^-)(^-|^ 12 |^-)( a3 (z)-P-)(z-Q*) 



(Ki-\$ M5 \Kt)(z-Q' 2 ) 



(Kl\t 123 \K b 2 -)(Kl\t 6 \K b 5 -)(a 3 (z) - P')(z - Q§) 



+ - 



Cll , CT6 (ifh# 456 |^-}(a 3 (z) - P 2 ') Cl2 , CT3 (^-|^ 456 |^-)(^-|^|^-)(a 3 (z) - P 2 *)(z - 



(^-|^ 6 |^-)( a3 (z)-P 3 -) 



(^-|^ 6 |^-)(KM^ 345 |^-)(a 3 (z) - P 3 ')(z - 



+ 



Cll ^(Kl\^ 12 \Kl-){z - Q' 2 ) c 12 , ff8 (^-|^ 6 |^-)(^-|^ 61 |^-)(a 3 (z) - Pf)(z - QJ 



(ifM^l^-)(^-o 3 ) 



(^-|^ 6 |^-)(KMfe 45 |^)(« 3 (z) - P 3 -)(z - Q5) 



(A.69) 



K^{z) = parity conjugate of K\{z) (obtained by applying eqs. ( 4. 12 )-( 4.17 1) 



(A.70) 



where 0:3(2;) can be read off from the on-shell values quoted below solution S\ in Fig. [3] 



K 3 {z) = c 7j0 



1 



c 23,cri 



1 



c 21,o-io 



4 {Kl-\h 23 \K\-){K\-\h\Kl-){z P 2 )(f3 3 (z) Q§) 2 (^-|^ 123 |^-)( z - P 2 ) 

1 c 23,<t 3 1 c 21,ct 3 

4 {K\-\UK\-){K{-\^ A5 \Kl-){z - P 3 ){(3 3 {z) - Q5) " 2 <^-|# 34g |i^-}(/3 3 (*) - Q») 



C21. 



1 c 24l(Tl (Ar I&12I *r>(&(*) - Q5) + ^<<*a >(* - Pa) 



2(^-|^ 6 |^-)(z-P 3 ) 2 (^-|^ 123 |^-)(^-|^|^)(z-P 2 )(/3 3 (z)-05) 



C21,, 



1 c 24 , CT7 (^-|^ 56 |^-)(z - Px) + c 24 , CTl0 (^-|^ 1 |^)(/3 3 (z) - Q\) 



2 (A^M^|^-)(/3 3 (z) - Q5) 2 (^-|^ 123 |^-)(xM^ 6 |^-)(z - P 2 )(/3 3 (z) - Ql) 

1 c 24 , ff3 (^-|fe 6 |^-)(z - PQ + c 24 , CT6 (^-|^ 61 |^-)(/? 3 (z) - Ql) 

2 (^-|^ 6 |^-)(^M^ 345 |^-)(z - P 3 )((3 3 (z) - Ql) 

1 c 24 , J8 (^-|^ 456 |^-)(z - P 2 ) + c 24 , CTll (^-|^ 6 |^-)(/? 3 (z) - Qf) 

2 (^-|^ 6 |^-)(^-|^ 345 |/^-)(z - P 3 )(ft(z) - Q' 2 ) 

c w . aw {Kt\h,\K\-){z - PQ c 13 , CTl (^-|fe 6 |A- 1 b -)(^-|^ 612 |^-)(z - P!)(/? 3 (z) - Qf) 
<K2"l*i23l*x">(* - ^) (^-|^ 123 |A^)(At|^ 6 |^-)(z - P 2 )(f3 3 (z) Ql) 

c w . a3 {Kl~\hi\K\')(Mz) ~ Ql) c 13>a3 {K\-\h 6 \K\-){Kl-\UKl-){z - Pi){P 3 {z) - Ql) 



(Kl-$ 3i5 \K b 5 -)((3 3 (z)-Ql) 



(K b 2 -\UK\~){K\-\t 3i5 \Ki-){z P 3 )(P 3 (z) - Ql) 
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c 10t<Te {Kl\h 6 \KlKz - Pi) c 13 ,„ 1Q {Kl\# 6 \Kl-){Ki-\$ 61 \K»-)(z - P 3 )(f3 3 (z) - Q\) 



(K b 2 -\h\K{-)(z-P 3 ) 



{Ki-\^\K\-)(KUUKt){z P2)(Mz) Ql) 



c 10 ,„ 7 (K b 4 -\hi\Kt)(f3 3 (z) - Q'i) C i3, CT6 (^-|^5 6 l^")(^"lfel^>(^ ~ Pi)<Jh{z) ~ Q'i) 



(KtitM-mw-QZ) 



(K b 2 -\h\Kl-)(Kl-\h, 5 \Kl)(z P 3 )(fo(z) - Q5) 



c lhr7l0 (K' 2 -\^\K b -)(z - P 3 ) c 12 ^K' 2 -\h 6 \Kn(Kl\hi\K~)^ ~ Pi)(fa(*) ~ Q'i) 



(Kl-\h 2 ,\K\-)(z-P 2 ) 



(Kl-\^\K{-){KUh\Kt){z P 2 )(P 3 (z) - Q' 3 ) 



cii^{Kt\^\Kt){k{z) - g|) ci 2taia (Kl$ 6 \Kl-)(Kl\t 612 \Kl)(z - P a )Q9 3 (z) - Qg) 
(xM^ 45 |^-)(/3 3 (z) - Q5) (J^-|jH ia3 |^-)(JrM^|^-)(z - P 2 )((3 3 (z) - Q«) 

cii, gg (^-|^ 45 6l^~)(^ - P*) ci2,<r s {Ki-\hm\K\-){Kl-\UKl-){z - P 2 )(p 3 (z) - Ql) 
{K\-\UK\-){z - P 3 ) (K\-\UK{-)(Kl\^\Ki-)(z - P 3 )(f3 3 (z) - Q-) 

cn, a7 {Kt\hi2\K\-){Mz) - Ql) c 12 , ffg (^-|fe 6 |XJ-}(^-|^ 61 |^-}( Z - Pi){(3 3 {z) - Ql) 



{Kl\UKl-){p 3 {z)-QD 



(xM^|x 1 b -)(xM^ 45 l^-)(z-p 3 )(/33W-Q^) 



(A.71) 



K±{z) = parity conjugate of K 3 {z) (obtained by applying eqs. (4.12 1-( 4. 17| |) (A. 72) 

where f3s(z) can be read off from the on-shell values quoted below solution 1S3 in Fig. [3) 



A.6.2 Heptacut #6 of the left hand side of eq. (4.1) 



The result of applying heptacut ^6 to the left hand side of eq. (4.1) is 

1 



i=i jTi 7=1 



(A.73) 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the four different kinematical solutions yields 



U A T e ( z ) 



for i = 1, ... ,4 



(A. 74) 
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A. 7 Heptacut #7 

Note that in this section we will leave the ± in 7^, P- , Qf etc. implicit and simply write 71, Pj, Qi 
etc. for notational simplicity. The heptacut considered here is defined by the on-shell constraints 
in eqs. (3.3)-(3.9) with the vertex momenta 



K A 



k\2 



K 2 
K 5 



hi 



K 3 
K 6 




. 



(A.75) 



Applying this heptacut to the right hand side of eq. (4.1 ) leaves the following linear combination 
of cut integrals 



Cie 



<7\ 



+ C 22,cri 




c 22 



0-9 




+ c 18,o-i 




CTg 




We define the spinor ratios 
Pi 



Pi 



(K\Kp 7i(ga+7l) 
2{K\K\) li-Silh 

(K\k{) 7l (5 2 +7i) 



p = [K b 2 k!] 5iS 2 (l+S , i/7i) 
5 2[i^fci] 7i 2 -5 , iS 2 



P 2 



Pfi 



IjM giS 2 (l+5i/7i) 
2[K\Kl\ lt-S 1 S 2 



P 3 



{K\K\\ giS 2 (l+gi/ 7 i) 



2\K\K\\ 



7?-5iS 2 



^g^^ fc --f^g 2 T^■ fc3 ^ (s " + " 3) -^ 

2(^-|Jf 3 |if{-) 



7l(S 2 +7l) K \, 



7f-SiS 2 



S l S 2 ( 5 l+7l) 
7l(7f-SlS 2 ) 



^•fc 3 -i(s 13 +S23)+V / A 



71 (S2 + 7i ) ^ ) [-gj ^ ] + gi S2 ( 1 +Si / 7i ) ^ ) [K\ K\] 
2(7i(S2+7i)(^^>[^JXil + SiS a (l+Si/7i)(^^S>[XS«S]) 



Q 2 



(A.76) 



and their parity conjugates 

[jW 71 (52+71 ) 



Pi 



P* 



2[K\K\\ 7i-5i52 

[K\k{\ 7i(5 2 +7i) 
2[K\k 1 ] 7i 2 -5i52 



p. = (K|fci) 5iS 2 (l+Si/7i) 
5 2(^fci) 7i 2 -5i52 



P 

-<2 



P 



P* 

J 7 



(g*K») 5i5 2 (l+5i/ 7 i) 
2<i^} 7i 2 -5i52 



5i5 2 (l+Si/7i) 
2(K\K\) 7i 2 -5i52 



7l(S 2 + 7l) 
7^~5lS 2 



^^ fc 3- !^! 2 2 ( ! l ^^ 2■ fc 3-|(sl3+. 2 3)-VA 



7l(7f~SlS 2 ) 



7l(S 2 +7l) i^b 



7f-SlS 2 



2{K{-\UK' 2 -) 

glg 2 (Sl+7l) 



7l(7f~SlS 2 ) 



^•fe 3 -i(si3+S 23 ) + VA 



2(K b -\£ 3 \K b -) 



Ql 



7i(52+7i)(^^)[^^]+5i52(l+5i/7i)(^^)[^^] n . 
2(7i(5 2 +7i)(^>[^^] + 5i5 2 (l+S 1 /7i)(^^>[^^]) ' ^ 



2(K\K\) 



where the discriminant appearing in Pq, P7, P' , P* is given by 

A f - 'v^- ~:: kX ■ fe, - SiS t + s s X ll] K 2 • fc 3 - 1(^13 + s 23 ) 



/ 71(52+71) I^D . jL 



7T 



- 4Sl5 #-"55 2 ( p • • **) ■ 



(A.77) 



(A.78) 
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This heptacut belongs to case III treated in Section 3.1, and there are thus four kinematical 
solutions (shown in Fig. Q. Parametrizing the loop momenta according to eqs. (3.12)-(3.13), the 



on-shell constraints in eqs. (|3.3|)-(|3.9|) are solved by setting the parameters equal to the values 

Pi =0 



_ 71(^2+71) 



"2 



■Si S2 (Si +71) 



1 



(A. 79) 



7i(SiS 2 -7?y ' 

and those given in Fig. jij The heptacut double box integral ii6,<ri is jEiELi/r^ 2 ^^) 



where 



' (32 72 ((5 a + 7l)(^-|^i|^-) + (l + I) {K\-\$\\K\-)) z(z - Q\)) 1 



for * = 1, 2 



32 72 ((5 2 + 7i)(^-|#i|i^-> + + <J^-||t* 2 |JfJ->) *(* - Qi 



for j = 3, 4 
(A.80) 



and where we recall that the ± in 7^ , ifj^, etc. have here been left implicit. 



A. 7.1 Heptacut #7 of the right hand side of eq. (4.1) 



The result of applying heptacut ^7 to the right hand side of eq. (4.1) is 

4 



SEE / dzJ?(z)Xt(z) 

± i=l JT i 



where the kernels evaluated on the four kinematical solutions are 



Kf(z) = Ci6, CTl 



1 



1 



C22,a 9 P* 



2 (^-|^ 3 |^-)(P- - P«)(P« - P«) 2 (^-|^|^-) ( p. - p.)(p. - p.) 

1 222,(7! 0:3(2) 1 022,^0:3(2) 



(A.81) 



(A.82) 



2 <^-|fe|i^-)Mz) - i^M*) - P 7 ') 2 (^-|^|^-)( a3 ( z ) - P;)(a 3 (z) P 5 ') 



c ls , ai (K\-\UKl)(a 3 (z) - P?)(a 3 (z) - P 2 *) c 18 , ff9 (^M# 6 |i^-)(a 3 W - P?)(a 3 (z) - PJ) 



<«*-|fe|^-)(a 3 («) - PS)(a 3 (z) P 7 ') 



{K\-\h\K\-){a,{z) Pl)(a,{z) P') 



K%{z) = parity conjugate of K^(z) (obtained by applying eqs. (4.12 1-( 4. 17 1) 
K&{z) — parity conjugate of K% {z) (obtained by applying eqs. (4.12 l-(4.17l) 
where a%{z) is given in Fig. El 



(A.83) 

(A.84) 
(A.85) 



A. 7. 2 Heptacut #7 of the left hand side of eq. (4.1) 



The result of applying heptacut ^7 to the left hand side of eq. (4.1) is 

4 



£W dzJt{z)\[A- 

± i=i JTi j=i 



(A.86) 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the four different kinematical solutions yields 

6 

= for 2 = 1,..., 4. (A.87) 



l[Aj™(z) 
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A.8 Heptacut #8 



This heptacut is defined by the on-shell constraints in eqs. (3.3)-(3.9) with the vertex momenta 

(A.88) 



Ki = ki 
K i = k 3 



K 2 = k 2 
K 5 = k 4 



K 3 = 



Applying this heptacut to the right hand side of eq. (4.1 ) leaves the following linear combination 
of cut integrals 




(h + hseYih + kse) 



+ C ll,cril — 




(h + kw) a 



+ C 24,cr 6 




+ c 24,ct 8 




C 2 4 



Til 



(li + kise) 
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We define the spinor ratios 

Kl-K 6 + \S (i +{Kl+K 6 yK\ 



Pi 
Pi 

Qi 



(Kl-\Jjtl+]jt 6 \Kl-) 

(K\-\$\+if: 6 \Kl-) 

2{K\K\)\K\K\\ 
K\-K (i + \S 6 + {K\+K & )-Kl 

(Kl~\t\+if; 6 \K b A -) 
and their parity conjugates 



P'i 



P* 



Ql 



(Kt\*\+*i\iqr) 
(kI-\$\+$ g \k\~) 

2\K\k\\{K\K\) 
K\-K 6 +\Ss+(K\+K 6 )-K\ 



■:\+$ 6 \kI-) 



Pi 



p* 

r 2 



Ql 



K\-k 6 



(J£*-[*b|K*-) 



(Kl-\h\KT) 



K\-k 6 



{K*~\h\K~) 



{K\-\h\Kl~) 



P 3 



p* 

^3 



Ql 



\k\-\^\k\~) 



^■fcs6+§S56 
<^-|^ 56 |4-> 



' (K b 4 -\H 56 \K b 5 -) 



(A.89) 



(A.90) 



This heptacut belongs to case I treated in Section 3.1 and there are thus six kinematical solutions 
(shown in Fig. [2|. Parametrizing the loop momenta according to eqs. (3.12 )-( 3. 13 ) , the on-shell 
constraints in eqs. (3.3)-(3.9) are solved by setting the parameters equal to the values 



«i = 1 
«2 = 



02 = 1 



(AM) 



and those given in Fig. [2] with 



(Ktm + telKtKz-Pi 



(A.92) 



2{K^Kl)[K\K\](z-P±) 
for kinematical solution 5s. The heptacut double box integral Ii7 >ai is Ylt=i §r ^ z ^i( z ) where 



Ji(z) 



1 



327172 



x < 



(Kl-\r 5 + te\Kl-)z(z-P^ 
(Kl-\t\ + t,\K\-)z{z-P l ) 
(K\-\t\ + $,\Kl)z{z-Q\) 
(Kl\$\ + t & \K\-)z{z-Q x ) 



for i = 2, 6 

for i = 4, 5 

for i = 1 

for i = 3 . 



(A.93) 



A. 8.1 Heptacut #8 of the right hand side of eq. (4.1) 



The result of applying heptacut ^8 to the right hand side of eq. (4.1) is 



1 6 / 

f dzJi(z)Ki(z) 

i=i 



(A.94) 



where the kernels evaluated on the six kinematical solutions are 



55 



Kl(z) = C 17 , CT1 + - t d.x 1 9 



C21,ct 3 



2 - P|) 2 (Xj-I^^-)^ - Q.) 

, c w , a3 (Kl$ 561 \Kl-)(z-Ql) 



c 23,<r 3 



4(^-|4|^)(^-|^|^->(P--P-)(z-Q-) (JfMfe|^-)( Z -Q«) 

1 c 2 ^ 3 {K\-\hs\K b 2 -){Pl - PI) + C24, g9 (^- fenfcKf - gj) + c 24 , CTll (^|fe 6 l^5')(^ - OS) 

2 (^-|^ 6 |^-)(^-|^|^-)(P- - p;)(z - Q5) 

ci 3 , (T 3(^-|^ 56 |^-)(^-|^ 56 |^-)(P 1 '-i^)(z-QS) Cll , g3 (^-|^ 56 |^-)(z-Qg) 
{K\-\UKl){Kl\h\Kl-){Pl-P^){z-Ql) {K\-\UKl-){z-QD 

c 12 ^(K\-% 6 \K^(K\-\h,i\Kl-){Pi - - Ql) , c n , CT11 <^-|^ 6 |^-}(P 1 - - P 3 *) 



(^-|^|^-)(^M^ 5 |^-)(P- - P-)(z - Q5) 



(JTM^|Jf*-)(i7-i^) 



(A.95) 



#2(2) = c i7,< 



1 



C21, 



1 



C21, 



1 



2 (XM^ 6 |^-)(Z - P-) 2 (^-|^ 5 |^-)(Q. - Q.) 



+ 



+ 



+ 



^{K\-\UKl-){K\-\UKl-){z-PS)(Ql-Ql) (K\-\^\Kl){z-PZ) 
I c^Kt^K^jz-P') + c 2i ^{K\-\UJKl-){z-Pl) + C24, gll (gI-^ 56 l4~)(Q?-93) 

2 (^M^ 6 |^-)(^M^ 5 |^-)(z - P 2 -)(QJ - Q«) 

C 1 3,a 3 (^-|fel^)^l^ 6 |^5"}(^-^3 , )(Q!-QI) , CU, CT3 (^-|^ 6 |^-}(QI-QS) 



(XM^ 5 |^-)(QJ-Q5) 



(^-|^ 6 |^-)(ifM^ 5 |^-)(z - P 2 *)(QJ - Q5) 



(^-|^ 6 |^-)(z-P-) 



(A.96) 



•^3(2) = parity conjugate of K\{z) (obtained by applying eqs. (4.12 l-( 4~T7| ) 



(A.97) 



Ki(z) — parity conjugate of K 2 {z) (obtained by applying eqs. (4.12 1-( 4.17| )) 



(A.98) 
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K 5 (z) = Cn jai 



c 23,a 3 



1 c 24 , CT3 (Kl \% m \K\-){z - P 3 ) + c 24 , ff6 |fe 61 I^5~>(& (*) 



(^-|^ 6 |^-)(ifM^ 5 |^)(z - P 2 )(/3 3 (z) - Q5) 



c 21,crii 



2(^-|^ 6 |^-)(z-P 2 ) 



l C2 4, CT8 (^-|^ 456 |^-)(z-F 1 ) + C2 4,. 11 (^|^ 56 |^-)(/33W-Q^) , 1 



C 21,<7 3 



2<^-|fe|^-)(^ 3 («)-Q5) 



(KM^ 6 |^-)(KMfe|^-)(z - P 2 )(/3 3 (z) - Q5) 



{Ki-\UK\-){z-P 2 ) 



ci 3 ^(^-|^ 456 |^-)(^-|fe 61 |^-)(z - Pi)Q8 3 (z) - g|) c 10 , CT3 (^-|fe 61 |^-)(/3 3 (z) - QJ) 
(^-|^|^-)(^-|^ 5 |^-)(z - P 2 )(/3 3 (z) - Q5) (Jfi-|jH B |^-)C8 8 («) - 05) 

c 12 , g3 (^-^ 456 |^-)(^-|^ 56 |^-)(z- P0(/3 3 ( Z ) - gf) | Cll ^(Kl$ 5e \Kl-)(p 3 (z) - Q%) 
{k\-\^\k\-){kI-\uk\-){z - P 2 )(k(z) - 05) (^-|^-)(/3 3 (z) - 05) 

c 12 , ff8 (xM^ 56 |^-)^-|^ 61 |^-)( Z -P 3 )(/3 3 ( 2 ) - QI) , c u , CT11 ^-|^ 56 |ifJ-)( 2 - P 3 ) 



(^-|^ 6 |^)(xM^ 5 |^-)(z - P 2 )(/3 3 (z) - 05) 



(^-|^ 6 |^-)(z-P 2 ) 



(A.99) 



Kq(z) — parity conjugate of K^{z) (obtained by applying eqs. (4.12)-(4.17)) 



(A.100) 



where /3^(z) is given in eq. (A.92). 



A. 8. 2 Heptacut #8 of the left hand side of eq. (4.1) 



The result of applying heptacut ^8 to the left hand side of eq. (4.1) is 



6 „ 6 

1=1 JVi 7 = 1 



(A.101) 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the six different kinematical solutions yields 



II 4 r 



3=1 



1 a tree 

16 



i 



i 



i 



Jl(z) \z z-Q? 2 



1 



X < 



Ja(z) \z Z-P-2 
J 6 (z) \z Z-P* 



for i = l 

for i = 4 

for i = 6 

for i = 2,3,5 



(A.102) 
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A.9 Heptacut #9 



This heptacut is defined by the on-shell constraints in eqs. ( 3.47 )-( 3.53 ). Applying it to the right 



hand side of eq. (4.1 ) leaves the following linear combination of cut integrals 




C23,, 






(h + h*)\h + k V2 f 



Cl3,< 




+ Ci2, CT1 



(il + fc 6 ) 2 (^2 + fel) 2 




(/l + fc 56 ) 3 (/2 + fcl) 



+ C12 



trio 




/c 3 £ 4 
(/ 1 + fc 6 ) 2 (/ 2 + fc 12 ) 2 



C 21, 




ClO,, 




C24,, 




+ cn 




C24,< 



(Il + fee) 2 




(I 2 + fci) 



We define the spinor ratios 



p _ [fc3fc2] 
ri - 2[fc 3 fcl] 



2 ( <fc2fc6)+ TO <A;2fc5) ) 
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^1 ~ 2(fc 4 fc 6 > 



P-7 
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(fcifcg) 

'2{k 2 k a ) 



ke-ki2 + ^s 12 

{kg 1^12^5 } 



Ph 



Q 3 



fcl-fc56 + f «56 
(&2 1^56 1^1 ) 



[fclfce] 
'2[&lk S ] 



(A.103) 
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and their parity conjugates 

p. _ {kjk2) 
r l - 2{k 3 k 1 ) ' 



^2 - 2 [fc 2 fc 6 ] 



[^4^5] 

2[fc 4 fc 6 ] 



Q5 



fc6'fcl2 + |^12 



r 3 



Ql 



(^1 1^561^2 ) 



(kike) 
'2<fcifc 5 > 



(A.104) 



This heptacut was treated in Section 3.2 and there are four kinematical solutions (shown in Fig. 



eqs. (3.47)-(3.53) are solved by setting the parameters equal to the values 



[6|. Parametrizing the loop momenta according to eqs. (3.54 )-(3.55 ), the on-shell constraints in 

(A.105) 



OL\ = 1 
«2 = 



/3i=0 

02 = 1 



and those given in Fig. [6] The Jacobian associated with the heptacut ( 3.47 )-( 3.53 ) is 



Hz) 



1 



1 



16si2S 4 5S 56 Z 



for i = 1, . 
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A. 9.1 Heptacut #9 of the right hand side of eq. (4.1) 



The result of applying heptacut #9 to the right hand side of eq. (|4.1|) is 

4 

1 



1 4 / 

- f dzJi{z)Ki{z 
i=i 



(A.107) 



where the kernels evaluated on the four kinematical solutions are 



Ki(z) 



2(k^\^ 3 \k 2 )(z-P^) 

(c 23 , CTl + 4c 13)CTl0 (ferine I ^2 ><*rl^xl*e>(* - p ') (Ql - Ql) 

+ ±c 12 , ai (kl\hs\kz)(K\h\K){z-PZ){Ql-Ql) 
+ 4c 12 , CTl0 (K\h\k2)(K\fi 12 \k 6 )(z - PS) (Ql - Ql) 
-2{kl\h\K){z - P?) (c 21 , CT7 + 2c w , a7 {k7 \h\K) (Ql - Ql) 

+ 2c 11 , ( , 7 <fe 5 -|^ 2 |fe 6 -}(Q 1 -Q^)) 

+ 2c 24 , ai (K\k , i2\K)(Q'i-Q') 

+ 2c 24 , CT4 (K\h\k^)(z-PS) + 2c 24 , ff7 (fcr |# 56 1 fc 2 " ) (z - F$ ) 

+ 2c Miai0 {k^\U^){Ql-Ql) 
+4c 13>CT1 <fcr|^ 6 |fc 2 -)<fc 5 -|^ 2 |fc 6 -)^ - p 3 *) (q; - Q3)) 
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K 2 (z) = - 



2(fcrifei*3-)(z-i7) 



1 



z-P 4 * 



2 V^ 1 "' Dl '^ ' ' v " i " u/ L """"' J (k 6 k 4 

x (c 23)CTl +4c 13>(Tl0 <fer|^ 6 |fc^)<fc6|^|A; 5 -}(2-P 2 , )(Qi-Q3) 
+ 4 Cl2 , CTl (fcr|^ 56 |fc 2 -)(fc 6 -|^ 1 |A ;5 -)(z - P 3 ') (Qj - Q a ) 
+ ^ 12 ^ B (k^\^)(k ( r\^ 12 \k^)(z - P') (Q x - Q a ) 
-2(fcrlfel^)(z - Pf) (c 21 , CT7 + 2c 10 , ff7 <^|# 1 |fcg ) (Qi - Q 3 ) 

+ 2 Cll , CT7 (V^ 12 |fc5 _ )(Qi-Q2)) 
+ 2c 24 , CTl (fc^|^ 12 |fc^)(Qi-Q 2 ) 

+ 2C 24 , CT4 (fcf Ifti l fc 2~) (* - P') + 2c 24 ,a 7 (fcf 1^56 1*2 > (* " 3?) 
+ 2024,^^1^1^) (Ql-Q 3 ) 

+ 4ci 3 , ffl fa" llBelAa ) |# 12 Ifeg) (z - PS ) (Qi - O2)) 



Ks(z) = parity conjugate of K 2 {z) (obtained by applying eqs. (4.12 1-( 4.17| ) 



Ki(z) = parity conjugate of K\{z) (obtained by applying eqs. (4.12 1-( 4.17 1) . 



(A.109) 

(A.110) 
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A.9.2 Heptacut #9 of the left hand side of eq. (4.1) 



The result of applying heptacut ^9 to the left hand side of eq. (4.1) is 

l 



:W dzJi(z)f[Af 

i=l JTi 1=1 



(A.112) 



where, assuming without loss of generality the external helicities are (1 ,2 , 3 + , 4 + , 5 + , 6 + ), the 
cut amplitude evaluated on the four different kinematical solutions yields 



i=i 



16 



-A 



tree 

h+H 



X < 



St 



1 / 1 



1 



J 2 (z) \z z-P> 
l/l 1 



Ja(z) \z z-P-2 





for i = 2 
for i = 4 
for £ = 1,3 
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